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ABSTRACT 

Motivated by the remarkably narrow range of measured spin frequencies of ~ 20 
accreting (and weakly magnetic) neutron stars in the Galaxy, Bildsten (1998a) conjec- 
tured that their spin-up had been halted by the emission of gravitational waves. If so, 
then the brightest persistent X-ray source on the sky, Scorpius X-l, should be detected 
by gravitational wave interferometers within ten years. Bildsten (1998a) pointed out 
that small nonaxisymmetric temperature variations in the accreted crust will lead to 
"wavy" electron capture layers, and the resulting horizontal density variations near 
e~ capture layers create a mass quadrupole moment. Neglecting the elastic response of 
the crust, Bildsten (1998a) estimated that even e~ capture layers in the thin outer crust 
can develop the quadrupole necessary to balance accretion torque with gravitational 
waves, Q22 ~ 10 37 - 10 38 g cm" 2 for accretion rates M ~ lCT 10 - 2 x 1CT 8 M Q yr" 1 . 

We present a full calculation of the crust's elastic adjustment to the density pertur- 
bations induced by the temperature-sensitive e~ capture reactions. We find that, due 
to the tendency of the denser material to sink rather than spread sideways, neglecting 
the elastic response of the crust overestimates, by a factor of 20 — 50, the Q22 that 
results from a wavy capture layer in the thin outer crust. However, we find that this 
basic picture, when applied to capture layers in the deep inner crust, can still generate 
Q22 in the necessary range, as long as there are <J 5% lateral temperature variations 

at densities in excess of 10 12 g cm -3 , and as long as the crustal breaking strain is 
high enough. By calculating the thermal flow throughout the core and the crust, we 
find that temperature gradients this large are easily maintained by asymmetric heat 
sources or lateral composition gradients in the crust. If the composition or heating 
asymmetries are independent of the accretion rate, then for M < 5 x 10~ 9 Mq yr -1 
the induced quadrupole moments have approximately the same scaling, oc M 1 ' 2 , as 
that necessary to balance the accretion torque at the same spin frequency for all M . 
Temperature gradients in the deep crust lead to a modulation in the thermal emission 
from the surface of the star that is correlated with Q 2 2- In addition, a ~ 0.5% lateral 
variation in the nuclear charge-to-mass ratio in the crust will also result in a Q22 
sufficient to halt spin-up from accretion even in the absence of a lateral temperature 
gradient. 

We also derive a general relation between the stresses and strains in the crust and 
the maximum quadrupole moment they can generate. We show under quite general 
conditions that maintaining a Q22 of the magnitude necessary to balance the accretion 
torque requires dimensionless strain a ~ 10~ 2 at near-Eddington accretion rates, of 
order the breaking strain of conventional materials. This leads us to speculate that 
accreting neutron stars reach the same equilibrium spin because they all are driven to 
the maximum Q22 that the crust can sustain. 

Key words: accretion - dense matter - radiation mechanisms: gravitational - stars: 
neutron - stars: rotation 
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1 INTRODUCTION 

Recent discoveries by the Rossi X-Ray Timing Explorer 
indicate that most of the rapidly accreting (M > 
10 _11 Af Q yr _1 ) and weakly magnetic (B < 10 11 G) neu- 
tron stars in our Galaxy are rotating in a narrow range 
of frequencies around v a ~ 300 Hz (van der Klis 1998). 
From both evolutionary considerations and their galactic 
distribution, we know that the neutron stars in these low- 
mass X-ray binaries (LMXBs) are relatively old systems (see 
Bhattacharya & Van Den Heuvel 1991) that have accreted 
enough angular momentum to reach rotation rates closer to 
the breakup frequency of £s 1.5 kHz. Hence, some mecha- 
nism must be found to halt the spin-up. One possible expla- 
nation is that the neutron stars have reached the magnetic 
spin equilibrium (where the spin frequency matches the Ke- 
plerian frequency at the magnetosphere) at nearly identi- 
cal spin frequencies. This requires that the neutron star's 
dipolar magnetic field strength correlates very well with M 
(White & Zhang 1997; Miller, Lamb, & Psaltis 1998). In 
an alternate scenario, Bildsten (1998a) suggested that these 
stars are rotating fast enough so that appreciable angular 
momentum is radiated away as gravitational waves (GW), 
allowing an equilibrium where the angular momentum added 
by accretion is lost to gravitational radiation. Equilibria via 
gravitational wave emission from rotational instabilities at 
much more rapid rotation rates had been postulated earlier 
by Papaloizou & Pringle (1978) and Wagoner (1984). 

The angular momentum loss rate from quadrupolar GW 
emission scales as v s , so that gravitational radiation ef- 
fectively provides a "wall" beyond which accretion can no 
longer spin the star up. Bildsten (1998a) estimated that 
"wavy" electron capture layers in the neutron star's crust, 
caused by a large-scale temperature asymmetry misaligned 
from the spin axis, can provide the quadrupole needed 
(Q22 ~ 2 x 10 38 g cm 2 at the Eddington accretion rate, 
MEdd = 2 x 10~ 8 Mq yr _1 ) to reach this limiting situation 
at v s ~ 300 Hz. Another possibility is GW emission from 
a continuously excited r-mode (i.e., a Rossby wave) in the 
neutron star core (Bildsten 1998a; Andersson et al. 1999), 
though Levin (1999) has shown that such steady-state r- 
mode can potentially be thermally unstable. More recently, 
Brown & Ushomirsky (2000) have shown that steady-state 
equilibrium between the accretion torque and r-mode grav- 
itational wave emission is incompatible with observations 
of the quiescent luminosities of neutron star transients. Fi- 
nally, Bildsten & Ushomirsky (2000) have shown that the 
extra dissipation due to a viscous boundary layer between 
the crust and the core is many orders of magnitude stronger 
that the viscous mechanisms previously considered, mak- 
ing it unlikely that the r-modes are excited in the cores of 
accreting neutron stars. We will not consider the r-mode 
hypothesis further here, but will instead concentrate on a 
self-consistent calculation of the mass quadrupole moment 
generated in the crust. 



1.1 Equilibrium Quadrupolar Gravitational Wave 
Emission 

We start by calculating the quadrupole moment that the 
neutron star (NS) must have so that the spin-up torque 
from accretion, N a , is balanced by emission of quadrupo- 



lar gravitational radiation. Consider a NS that is per- 
turbed from sphericity by a density perturbation Sp = 
Re{8 pim(r)Yim(d , 4>)}- Let Qi m be the perturbation's multi- 
pole moment, defined by 



Qi 



5pim(r)r l+2 dr. 



(1) 



We concentrate on perturbations Q22 with I — m — 2; these 
perturbations radiate at a frequency u gw = 2v s . The result- 
ing rate of loss of angular momentum N gw is 



_ E gw _ G 1 ,-■ --ab 
3W ~ O _ r-5 ZO\ labl I 



(2) 



Q c 5 50. 

where I ab = J Spr a r b d 3 V, Q = 2tiv s , and . \" means 
"time-averaged over one period." A little algebra shows that 

(lab I ) = (2567r/15)r2 6 Q22! so equation (2) becomes 
_ 2567T GtfQh 

For simplicity, we assume that the accreted angular mo- 
mentum is that of particles arriving from the inner edge 
of the accretion disk (placed at the NS radius), so that 
N a = M(GMR) 1/2 , where M and R are the mass and ra- 
dius of the NS. The required quadrupole moment Q cq such 
that gravitational wave emission is in equilibrium with the 
accretion torque, N a , is then 

Q cq = 3.5 x 10 37 g cmMl^R\ /A 

M V /2 ( 300 Hz" 



(4) 



1O" 9 M yr- 



5/2 



where M 1A = M/1AM Q and R 6 = R/W km. The range of 
M's typically encountered in the low-mass X-ray binaries is 
10~ 10 -2x 10~ 8 M Q yr-\ requiring Q22 ~ 10 37 -10 38 g cm 2 
for v s = 300 Hz (Bildsten 1998a). 

This paper is devoted to learning whether a quadrupole 
this large can be generated from deformed capture layers 
in the crust (Bildsten 1998a), and to finding the magnitude 
and distribution of the elastic strain required to sustain it. 
We answer the latter question in some detail in § 6, where 
we show that Q22 is related to the typical strain, a, via 

Q 2 -i- 2 xio 3 w(^)|g. 



(5) 



A more complete version of this relation that includes the 
dependences on the density at the crust-core transition and 
the composition of the crust is given in equation (69) . If we 
presume that the crust is pushed to some yielding strain 
o"max by the physical effects calculated here, we find 



295 Hz 



10" 



2/5 M°d 



M 



10" 8 M yr- 



1/5 



, (6) 



a rather suggestive relation that points to a possible answer 
as to why so many LMXBs are in a narrow spin frequency 
range, namely that accretion always drives the crustal strain 
to the breaking point. 

Before we launch into the detailed discussion of how 
the crust is stressed and how it responds, we calculate the 
strength of the GW signal from such an equilibrium radia- 
tor. Consider a neutron star at a distance d with an energy 
flux in gravitational waves, E gw , but with an unknown spin 
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orientation to the observer. We define the source's "angle- 
averaged" field strength h a (at Earth) by 



h l = ^j (h 2 + + hl)dQ 



(7) 



where the integral is over source orientations . The standard 
formula for the effective stress-energy of gravitational waves 
yields 



cWfo 2 
G 



(8) 



which, when we write it in terms of Q22 (using equation 3) 
gives 

16 /7r\V 2 GQ 22 fi 2 



16 /7TV 

K = T UJ 



dc 4 



(9) 



When the angular momentum loss by gravitational radiation 
balances angular momentum gain by accretion, 



^3, X 10-^« 1/2 
Ml'* V v s ) 

1/2 



(10) 



10 -8 erg cm 2 s 



(Wagoner 1984; Bildsten 1998a). Here we have replaced M 
and d with the observed X-ray flux, F x = GMM/kirRd 2 . 
The gravitational wave strength for a neutron star accret- 
ing at the Eddington limit at the galactic center is then 



h a 



5 x 10" 



Prior accurate knowledge of the position 



on the sky and orbital periods of many of these X-ray bi- 
naries will allow for deep searches with the suite of laser- 
interferometric gravitational wave detectors currently under 
construction (LIGO, VIRGO, GEO-600, and TAMA-300; 
see Brady et al. 1998 and Brady & Creighton 2000). The 
nearby source Scorpius X-l is the obvious first target. Its 
X-ray flux is F x « 2 x 10~ 7 erg cm~ 2 s~ x and the spin pe- 
riod is still unknown, but likely in the range of 300 Hz (van 
der Klis 1998) giving h a « 1.6 x 10~ 26 . Bradshaw et al. 
(1999) recently determined that Sco X-l is at a distance 
d = 2.8 ± 0.3 kpc, giving a luminosity close to the Ed- 
dington value for cosmic abundances, 2 x 10 38 erg s _1 , im- 
plying M ~ 2 x 10~ 8 Mq yr" 1 and a required quadrupole 
Q22 « 1.6 x 10 38 g cm 2 (300 Hz/z/ s ) 5/2 for equilibrium GW 
emission. 



1.2 Origin of the Crustal Quadrupole Moment 

The crust of a neutron star is a thin (« 1 km) layer of crys- 
talline "ordinary" (albeit neutron-rich) matter that overlies 
the liquid core composed of neutrons, protons, and electrons 
(see Figure 1). The composition of the crust, i.e., the mass 
number A and the charge Z of the nuclei, varies with depth. 
As an accreted nucleus gets buried under an increasingly 
thick layer of more recently accreted material, it under- 
goes a series of nuclear reactions, including electron cap- 
tures, neutron emissions, and pycnonuclear reactions (Sato 
1979; Haensel & Zdunik 1990b; Blaes et al. 1990). The crust 



* A common measure of GW source strength is the characteristic 
amplitude h c . It is related to h a by h c ~ l.l5h a . 



Accretion 




Figure 1. The schematic structure of the neutron star, displaying 
the inner and outer crust lying on the liquid core. The ocean 
on top of the crust ends when the material is dense enough to 
crystallize. The outer crust consists of nuclei in a lattice and a 
background of degenerate, relativistic electrons, whereas the inner 
crust (at densities > (4 — 6) x 10 11 g cm -3 ) has free neutrons in 
addition to the nuclear lattice. 



therefore consists of layers of different nuclear composition, 
as indicated schematically in Figure 1. 

Essentially all the pressure in the outer crust, and an 
appreciable fraction of the pressure in the inner crust, is sup- 
plied by degenerate relativistic electrons. The pressure must 
be continuous across the boundaries of the compositional 
layers. However, electron capture reactions reduce the num- 
ber of electrons per nucleon, and hence require density jumps 
between the compositional layers depicted in Figure 1. In the 
outer crust, these density jumps are as large as « 10%, while 
in the inner crust the density contrast is smaller, < 1%. 

Bildsten & Cumming (1998) showed that at the typical 
crustal temperatures of accreting neutron stars, >2x 10 8 K, 
the electron capture rates are sensitive to the local tem- 
perature (see discussion in § 2.1). In this case, regions of 
the crust that are hotter undergo electron-capture transi- 
tions at a lower density (and thus larger radius) than the 
colder regions. This effect is illustrated schematically in Fig- 
ure 2, which shows a patch of the crust near a capture layer 
where nucleus (Ai, Z\) is transformed into (A2, Z2). If there 
is no lateral temperature gradient, then the capture layer is 
spherically symmetric, as indicated by the dashed line. The 
lateral temperature gradient (shown by the arrow) induces 
some extra captures on the left (hotter) side, and reverses 
some captures on the right (cooler) side. If the lateral tem- 
perature contrast is ST, then e~ captures on the hot side of 
the star happen at the Fermi energy that is roughly YksST 
lower, and captures on the cold side proceed at the Fermi 
energy about TksST higher than in the unperturbed case 
(Bildsten 1998a). Typical values of T are 10 — 20 (see the 
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Ocean 




Figure 2. A cartoon description of how a transverse temperature 
gradient in the crust will lead to a varying altitude for the electron 
captures. The dashed line denotes the unperturbed location of the 
e — capture boundary between the layer of material composed of 
nuclei (Ai,Zi), a mass fraction X n i of free neutrons, and free 
electrons, and a layer with nuclei (Az,Z2), X n 2 fraction of free 
neutrons, and free electrons. When a lateral temperature gradient 
VT is introduced, the capture boundary shifts to a new location 
shown by the solid line. Note that this cartoon presumes that the 
crust is infinitely rigid and therefore does not adjust elastically 
to the shift Azd of the capture boundary. 



top panel in Figure 3) . Since the Fermi energy increases with 
depth, captures occur a distance Azd higher in the crust on 
the hot side, and a distance Azd lower on the cold side. 
The resulting capture layer is indicated by the solid line in 
Figure 2. 

Therefore, a large-scale temperature asymmetry de- 
forms the capture layers. Such a temperature gradient, if 
misaligned from the spin axis, will give rise to a nonaxisym- 
metric density variation and a nonzero quadrupole moment 
Q22 even if the composition of the crust has no lateral vari- 
ation. Of course, a lateral composition gradient in the NS 
crust can exist without a temperature gradient, and would 
also create a quadrupole moment. 

On dimensional grounds, the quadrupole moment gen- 
erated by a temperature-dependent capture boundary is 
Q22 ~ Qfld = ApAzd^? 4 , where Ap is the density jump 
at the electron capture interface. This fiducial value, Qfld, is 
the quadrupole moment that would result if crustal matter 
just moved horizontally to regain horizontal pressure bal- 
ance. Using this estimate, Bildsten (1998a) argued that a 
single wavy capture boundary in the outer crust could gen- 
erate Q22 sufficient to buffer the spinup due to accretion, 
provided that temperature variations of ~ 20% are present 
in the crust. 

However an important piece of physics is missing from 
this fiducial estimate: the shear modulus of the crust This 
must be important to the estimate, since if \i vanishes, the 
crust becomes a liquid and cannot support a non-zero Q22- 
In addition, a NS crust is thought to have a small shear 
modulus relative to the pressure, n/p ~ 1(T 2 - 10" 3 typi- 
cally. Much of this paper is concerned with improving upon 
the above fiducial estimate. Treating the crust as an elastic 



solid, we assume the existence of a horizontal temperature 
or composition gradient, and then solve for the displacement 
field that brings the crust into equilibrium, with the gravita- 
tional, pressure, and shear-stress forces all in balance. From 
the density perturbation Sp we calculate Q22 for a single 
electron capture layer and find that it is typically 5 — 50 
times smaller than the fiducial estimate, depending on the 
capture layer depth (see § 5.1). Hence, in order to generate 
a Q22 large enough for substantial GW emission, we must 
include capture layers in the more massive neutron-rich part 
of the crust. 

Before launching into detailed calulations, we now de- 
scribe why it is plausible to assume that lateral temperature 
and composition gradients are present in the crusts of ac- 
creting neutron star. We present a detailed model of thermal 
gradients in the crust in § 3. 

1.3 Possible Causes of Temperature Asymmetries 

In LMXBs, accretion will replace the primordial crust after 
about ~ 5 x 10 7 yr (M/10~ 9 M© yr" 1 ). This is much shorter 
than the lifetime of such systems. While it is plausible that 
the primordial crust is spherically symmetric in composition, 
we suggest that the accreted crust need not be. 

The accreted crust is composed of the compressed prod- 
ucts of nuclear burning of the accreted hydrogen and helium 
in the NS's upper atmosphere. The nuclear mix entering the 
top of the crust depends sensitively on the burning condi- 
tions and is still not well known. Schatz et al. (1999) showed 
that the products of steady-state burning in the upper atmo- 
sphere are a complicated mix of elements far beyond the iron 
peak. The exact composition (and the average A) depends 
on the local accretion rate, which could have a significant 
non-axisymmetric piece in the presence of a weak magnetic 
field. 

However, except in the highest accretion rate LMXBs, 
nearly all of the nuclear burning occurs in Type I X-ray 
bursts, sudden consumption of fuel accumulated for hours 
to days prior to ignition. The rotational modulation ob- 
served during type I X-ray bursts in several LMXBs (first 
detected in 4U 1728-34 by Strohmayer et al. 1996; see van 
der Klis 1999 for a review) provides conclusive evidence that 
bursts themselves are not axisymmetric. Until the origin of 
this symmetry breaking is clearly understood, it is plausi- 
ble to postulate that these burst asymmetries get imprinted 
into the crustal composition or result from them. Finally, it 
is possible that there is a feedback mechanism that causes 
composition asymmetries to grow: lateral composition vari- 
ations lead to temperature asymmetries (as shown below), 
which, in turn, affect the burning conditions of the elements 
entering the crust. 

Clearly, it is almost impossible to compute the mag- 
nitude of the composition asymmetry from first principles, 
and we shall not attempt to do so here. Instead, we postu- 
late that such asymmetry exists at some level and explore 
its consequences. In particular, we show that a composi- 
tion asymmetry in the crust will modulate the heat flux 
through it. Thus, one of our predictions is a relation between 
the modulation in the persistent thermal emission of accret- 
ing neutron stars and their quadrupole moments. In other 
words, the lateral temperature gradient in the deep crust 
that generates a Q22 sufficient to halt accretional spin-up 
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also gives rise to a certain thermal flux modulation, which 
we quantify in § 3.3 and § 5.2. 

The consequences of the lateral composition asymme- 
try of the crust are two-fold. First, different elements have 
different charge-to-mass ratio, and hence different thermal 
conductivity and neutrino emissivity, both of which scale 
as Z 2 1 A. This lateral variation of the transport properties 
modulates the heat flux in the NS crust, leading to lateral 
temperature variations ST. Secondly, the nuclear reactions 
in the deep crust release energy and heat it locally. The 
heat release is again dependent on the particular element, 
and hence varies laterally if the crust is compositionally 
asymmetric, also giving rise to a temperature gradient. We 
quantify both of these thermal effects in § 3. In addition 
to these thermal effects, a composition gradient generates a 
quadrupole moment directly, because elements with differ- 
ent Z/A have different characteristic electron pressures, and 
the resulting transverse pressure gradient elastically deforms 
the crust. 



1.4 Outline of the Paper 

The purpose of this paper is three-fold. In the first part, 
we calculate the temperature asymmetry in the crust that 
can arise from asymmetric conductivities and/or heating. 
As stressed in § 1.3, at this point it is just a (testable!) con- 
jecture that there are large-scale temperature asymmetries 
misaligned with the spin of the neutron star. We show how 
these can arise due to composition gradients, despite the 
large thermal conductivity of the core. In § 2, we discuss 
the structure of the neutron star crust; in particular, we ex- 
tend the work of Bildsten (1998a) and Bildsten & Cumming 
(1998) on the structure of electron capture layers by consid- 
ering regions where there is a large density of free neutrons. 
In § 3, we calculate the magnitude of the temperature vari- 
ations induced by the composition variations, both through 
their effect on the conductivity and on the local heating 
rate in the crust. We find that 10% lateral variations in the 
heating rate or conductivity result in < 5% temperature 
asymmetry in the deep crust at accretion rates of order the 
Eddington rate, and <; 1% lateral temperature variations at 

KT 2 M E dd. 

In the second part of our paper (§ 4 and 5), we calculate 
the elastic adjustment of the crust induced by a deformed 
electron capture layer or a smooth composition gradient, and 
determine the resulting mass quadrupole moment, Q%2- If 
we consider the Q22 generated by a single deformed capture 
layer in the outer crust (which contains only ~ 1O~ 5 M ), 
we find that Q22 is smaller than the fiducial estimate of 
Bildsten (1998a) by a factor 20 — 50. But deformed capture 
boundaries in the deep, inner crust can generate sufficient 
Q22 to halt the spin-up of neutron stars at 300 Hz, provided 
there are ~ 1% lateral temperature variations, and provided 
the induced strains do not crack the crust. Quadrupole mo- 
ments due to multiple capture layers add linearly, and hence 
the required temperature asymmetry is even smaller. More- 
over, a smooth 0.5% composition gradient results in a sim- 
ilar quadrupole moment even in the absence of a lateral 
temperature gradient. Our solutions exhibit typical strains 
a > 10~ 2 at near-Eddington accretion rates, and a ~ 10~ 3 
at 10~ 2 MEdd (see § 6.3). The former is larger than the break- 



ing strain for terrestrial rocks under atmospheric pressure, 
but is perhaps possible for highly compressed neutron star 
crusts. The level of crustal strain required to sustain the 
large quadrupole moments is perhaps the most problematic 
feature of our model. 

The third part of this paper is an investigation of the 
relation between Q22 and crustal shear stresses. Specifically, 
we derive a relation, Eq. (64), that expresses Q22 as an in- 
tegral over shear stress terms in the crust. Eq. (64) holds 
independent of any detailed model of how those stresses are 
generated, and it immediately gives us an upper limit on Q22 
for a given crustal breaking strain. In § 6.5, we estimate the 
correction to the results due to the neglect of the gravita- 
tional potential perturbation (the Cowling approximation). 

Finally, we close in § 7 with a summary of our efforts 
and a discussion of what is still missing from the theoretical 
picture. 



2 STRUCTURE OF THE ACCRETED 
NEUTRON STAR CRUST 

Bildsten (1998a) confined his discussion to the outer crust 
(before neutron drip at p < p n d ~ (4 — 6) x 10 11 g cm~ 3 , 
Haensel & Zdunik 1990b), which is held up by relativistic de- 
generate electrons. However, as we found in the course of our 
work, the capture layers that produce the largest quadrupole 
moments are at densities much greater than neutron drip, 
and so we need to model the entire crust. Our modeling of 
the thermal structure of the crust and core mostly follows 
Brown (2000). 

The nuclear mix entering the top of the crust is not 
well known. This mix depends sensitively on the conditions 
of hydrogen and helium burning in the upper atmosphere. 
For steady burning, Schatz et al. (1999) showed that the 
products are a complicated mix of elements far beyond the 
iron peak. The more relevant case of time-dependent nu- 
cleosynthesis in X-ray bursts is still unresolved (Rembges 
et al. 1997; Thielemann et al. 1998; Rembges et al. 1998). 
What eventually needs to be done is a calculation of the nu- 
clear evolution of these complicated mixes throughout the 
deep crust. Lacking these inputs, we take for our model 
of the crustal composition the tabulation given by Haensel 
& Zdunik (1990b, a). These authors start by assuming that 
pure iron enters the top of the crust, and that at each pres- 
sure the crust is composed of a single nuclear species, with 
the transitions between species being abrupt. Starting with 
56 Fe, their calculations produce the sequence of the most en- 
ergetically favorable nuclei (and the range of densities and 
pressures for each one) under the constraint that only elec- 
tron captures, neutron emissions, and pycnonuclear reac- 
tions are allowed (i.e. the crust is too cold for a thermal 
reshuffling of nucleons to occur). 

Taking the composition to be pure instead of mixed 
has negligible effect on the hydrostatic composition of the 
crust. However the composition does have a large effect on 
the crust's thermal conductivity, K; we use an estimate of 
K that is appropriate for a lattice of mixed composition 
(Schatz et al. 1999). Also, while Haensel & Zdunik (1990b,a) 
approximate the capture transitions as infinitely sharp, we 
resolve the actual, finite-thickness capture layers by integrat- 
ing the capture rate equation, following Bildsten & Cum- 
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ming (1998) . We extend their work by including the presence 
of free neutrons, which allows us to resolve capture layers at 

P > Pnd- 



2.1 Temperature Sensitivity of Electron Capture 
Rates 

In accreting neutron stars, the nuclear transformations of 
the crustal material are driven by increasing compression 
from the weight of the overlying matter. The location and 
thickness of the reaction layers are determined by the com- 
petition between the corresponding reaction rate (at a given 
p and T) and the local compression timescale, i.e., the rate 
at which the local conditions are changing. The compres- 
sion timescale at a given depth is f comp = p/rhg, where 
rh = M/4irr 2 is the local accretion rate, g = GM r /r 2 is the 
local gravitational acceleration, and p/g is approximately 
the column density, J r pdr. This is just the time it takes 
for the pressure on a fluid element to double due to the extra 
hydrostatic pressure of new material added at the top. 

Consider the transformation of a region of the crust 
where the predominant nucleus has charge Ze and mass 
Am^. (For simplicity, consider a region where p < p n d-) 
As it is compressed, the electron Fermi energy, Ef, rises 
to the point where an electron capture on the nucleus is 
energetically allowed. The reaction transforms an element 
(A, Z) into (A, Z — 1). In practice, the mass difference be- 
tween (A, Z — 1) and (A, Z — 2) is always greater than that 
between (A,Z) and (A, Z — 1), so the subsequent reaction 
(A, Z — 1) + eT — > (A, Z — 2) + u e is very fast, and typi- 
cally proceeds immediately (Haensel & Zdunik 1990b; Blaes 
ct al. 1990). The second capture releases of order 1 MeV, 
making the process effectively irreversible. The two succes- 
sive electron captures can often be treated as one reaction, 
with the rate-limiting step being the first capture. At densi- 
ties greater than neutron drip, the captures are accompanied 
by neutron emission (typically 6 neutrons are emitted), and 
sometimes also by a pycnonuclear reaction. 

What is the exact place where the reaction becomes fast 
enough to compete with compression? At T = 0, in order 
for electron captures to proceed, the fluid element must be 
compressed until the electron Fermi energy Ef exceeds the 
threshold energy Q (roughly the mass difference between 
the e~ capturer and the product, Blaes et al. 1990). Hence, 
at low T the electron capture rate is not very temperature- 
sensitive, and the location of the capture layers does not 
depend on the local temperature. The thickness of the elec- 
tron capture layers in this case is set by the need to have 
enough phase space (determined by Ef — Q) so that the 
electron captures proceed at a rate comparable to that of 
the compression (Blaes et al. 1990). 

However, Bildsten & Cumming (1998) showed that 
when the temperature is in excess of 2 x 10 8 K (condi- 
tions typical for crusts of neutron stars in LMXBs), there are 
enough electrons on the thermal tail of the Fermi-Dirac dis- 
tribution so that captures can proceed even when Ef < Q- 
The capture rate R ec is approximately given by 



Re 



In 2 
It 



2Q 2 (k B Tf 
(m e c 2 ) 5 



■ exp 



V k B T ) 



(11) 
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Figure 3. Top: Sensitivity T = —AEF/k B 8T of the location 
of c~ capture layers to local temperature perturbations as a 
function of the capture layer threshold Q. Bottom: Vertical shift 
of a capture layer Az^ for a fiducial temperature perturbation 
ST/T = 0.05, disallowing elastic readjustment of the crust, plot- 
ted as function of the capture layer threshold energy Q. These 
plots are generated using the formalism described in Appendix A 
(see Eqs. [A18] and [A20]). The complicated dependence of Az^ 
and T on the capture layer depth stems from the differences in 
the number of electrons captured and number of neutrons emitted 
in each capture layer. 



ft value is for the first e~ capture transition and depends 
on the degree of forbiddenness of the reaction. Typical ft 
values range anywhere from 10 3 to 10 s seconds for the tran- 
sitions that are relevant. In this paper we use ft = 10 4 sec 
for all reactions. As we show below, the dependence of the 
location of the capture layer on ft is only logarithmic, so 
even an error of a few orders of magnitude is not important 
to our work. On the other hand, the sensitivity to the local 
temperature is exponential, which is why even modest lat- 
eral temperature gradients can generate sizeable quadrupole 
moments. 

How does the location of the capture layer change 
with the local temperature? Most electron captures hap- 
pen when the lifetime of an element to electron capture, 
t cc = 1/Rcc becomes comparable to the local compression 
timescale t comp (Bildsten & Cumming 1998; Bildsten 1998a). 
Prethreshold captures then proceed at Ef ~ Q — Tk B T 
(Bildsten 1998a), where 



T = In 



p c (E F )\n2 2Q 2 (k B T) 3 



rag ft (m e c 2 ) 6 



10-20, 



(12) 



(equations [2] and [5] of Bildsten & Cumming 1998). The depending on the capture layer pressure p c and other param- 
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eters. The procedure for solving this equation is described 
in Appendix A. The function T is plotted in the top panel 
of Figure 3. 

We define Az d to be the distance that the capture layer 
would shift vertically if the crust were absolutely rigid; i.e., 
if there were no elastic readjustment. Then Az d is approxi- 
mately given by 



A£d 

h 



T 



dlnp k B T ST 
dlriEr^^' 



(13) 



where h is the local scale height. At densities lower than 
neutron drip, p oc E F , while for p > p nt j the dependence 
of p on Ef is even steeper, as electrons supply only a frac- 
tion of the total pressure. This relation is formulated more 
precisely in Appendix A, and Azd for a fiducial 5% temper- 
ature perturbation 8T/T is shown in the bottom panel of 
Figure 3. 



2.2 Capture Layers at Densities Higher than 
Neutron Drip 

An exact calculation of the crust's composition would re- 
quire a reaction network large enough to allow for the pos- 
sibility of several elements undergoing captures at the same 
time. This is beyond the scope of our paper. For simplic- 
ity we assume that between the capture layers, the crust is 
composed of a single species of nucleus (A,Z), electrons, 
and neutrons (once p > p n d)- Capture layers, however, 
contain a mix of elements: a mass fraction Xi of element 
(A\,Z\) that is abundant above the capture layer (see Fig- 
ure 2), a mass fraction X n of free neutrons, a mass frac- 
tion X2 — 1 — X\ — X n of clement (A2, Z2) (the end prod- 
uct of electron captures onto (Ai, Zi)), and electron density 
n e = p/ p c m,h. With this simplification, we only need to in- 
tegrate one rate equation at a time. 

We consider a capture layer, where elements (Ai,Zi) 
are transformed into (^2,^2), and, at the top of the layer, 
where the mass fraction of (A2, Z2) is zero, the neutron mass 
fraction is X n \. Suppose that the reaction is a simple one, 
i.e. it consists of capturing Z\ ~ Z2 electrons, and emitting 
A\ — A2 neutrons. Then, at the bottom of the layer, where 
the reaction is complete and the mass fraction of (Ai,Zi) 
is zero, the mass fraction of free neutrons X n 2 is such that 
j4i/(1 — X n i) = ^2/(1 — X n 2)- Similarly, if electron captures 
and neutron emissions in the layer are accompanied by a 
pycnonuclear reaction (i.e fusing of two nuclei (A, Z) into 
a single (2A,2Z) nucleus), then simple bookkeeping shows 
that - A nl ) = (1/2)A 2 /(1 - X n2 ). 

Now consider some point in the capture layer where 
both reactants (Ai,Zi) and products (^2,^2) are present. 
In practice, neutron emissions are triggered by electron cap- 
tures, and the timescale for neutron emission is always 
much shorter than the electron capture timescale (Sato 1979; 
Haensel & Zdunik 1990b). We therefore assume that these 
steps proceed simultaneously, so the proportion of neutrons 
per nucleus of each type stays fixed; e.g., if only {Z\ — Z2)/2 
electrons have been captured, then exactly (Ai — A2)/2 neu- 
trons have been emitted. This gives 



X n = 



AlAnl + (1 — A n l — Xl)X n 



i-x nl 

The electron mean molecular weight is 



(14) 
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Figure 4. Structure of an electron capture layer in the outer 
crust, at threshold energy Q = 23 MeV. The reaction is 
56 Ca+2e~ ^ 5e Ar+2u e . The upper panel shows the run of den- 
sity with pressure, with the density jump Ap/p fs 10% clearly 
evident. The middle panel shows the abundance X\ of 56 Ca. The 
bottom panel shows the electron Fermi energy Ep. Note that the 
slope of _Bp does not change in the capture region and is always 
dlnE F /dlnp = 1/4. 



(15) 



_L - XlZl I (l-^n-Al)Z 2 
Pc Al A2 

The change in mass fraction X\ is computed from the con- 
tinuity equation for species (Ai,Zi), 



+ V ■ (nXiv) = -nX 1 R cc {p,X 1 ,T), 



(16) 



d{nX 1 
dt 

where n is the baryon number density, 7i cc is the electron 
capture rate (11), and v = —fM/4irr 2 p is the very slow 
downward motion of the fluid clement. We assume that the 
accretion flow is predominantly radial and time independent, 
so Eq. (16) becomes 



dlnXi 
dr 



pRecfaXuT), 



(17) 



making clear the competition between the capture rate and 
the compression timescale, t cornp , defined earlier. 

We use the sequence of reactions shown in Tables 1 
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Figure 5. Top three panels are same as Fig. 4, but for a capture 
layer in the inner crust, with threshold Q = 95 MeV. For com- 
parison, dotted lines show what the run of p (top panel) and Ef 
(third panel) would be if the capture layer was not present. The 
fourth (bottom) panel shows the slope dlnSp/dlnp, computed 
using Eq. (A17). 

and 2 of Haensel & Zdunik (1990b) to describe the crust's 
composition. Their sequence starts with 56 Fe with X n = 
for densities p < 1.5 x 10 9 gr cm -3 and ends with 88 Ti with 
X n = 0.8 at densities p > 1.61 x 10 13 gr cm" 3 . For the 
threshold energy Q, we use the value of the electron chemical 
potential at the (abrupt) transitions between nuclei in the 
calculation of Haensel & Zdunik (1990a), as shown in their 
Table 2. Haensel & Zdunik (1990b) stop their calculations 
at p > 10 13 gr cm -3 , where their EOS state for accreted 
matter approaches the standard Baym et al. (1971) equation 
of state for cold catalyzed matter, and so the authors stop 
their tabulation. At these densities the equation of state is 
dominated by neutron pressure, and the exact A and Z of 
the nuclei have little effect on the EOS. But an important 
difference emphasized by Haensel (1997) is that Z f> 50 
at the bottom of a cold catalyzed crust, while the accreted 
model has Z ~ 20. While this has little effect on the EOS, 
the difference in Z does affect the shear modulus (Sato 1979; 
Haensel 1997). 

The crust extends to densities p « 2 x 10 14 gem -3 , 



and there are many more capture layers in the deep crust, 
each of which contributes to the NS quadrupole moment. 
Because the tables in Haensel & Zdunik (1990b) only extend 
to p — 1.61 x 10 13 gr cm -3 , in all our calculations we insert 
an extra, ad hoc, movable capture layer in the bottom part 
of the crust, and study the quadrupole moment induced by 
this layer as function of its position. We take the crust to be 
made up of (A, Z) = (88, 22) above the capture layer and 
(A, Z) — (82, 20) below. We selected these values because 
in most of the capture reactions listed by Haensel & Zdunik 
(1990b), two electrons are captured and six neutrons are 
released. 

Figures 4 and 5 show the structure of the capture lay- 
ers at densities below and above p nt j. The top panel shows 
the run of density with pressure (downward direction is to 
the right). At p < p nd (Figure 4), the pressure is supplied 
entirely by degenerate electrons. In this particular capture 
layer, Z changes from 20 to 18, so the density jump is 
Ap/p = Afi e /p e « 2/Z = 10%. The density change in a 
capture layer at p > p nd (Figure 5) is much smaller, as elec- 
trons provide a much smaller fraction of the pressure. 

For p > p nd , the capture layers become much thicker, 
both in pressure coordinates and in physical coordinates. We 
understand this as follows. Since the reaction rate (11) is ex- 
ponentially sensitive to Ef—Q, the transition is always sharp 
in Ef coordinates (see bottom panel of Figures 4 and 5). 
However, in pressure coordinates, the width of the capture 
layer is set by Ap/p = (din Ef / dlnp) -1 (AEf / Ef) , where 
cLIuEf /dlnp is given by equation (A17) in Appendix A. 
Exponential sensitivity of the reaction rate ensures that 
AEf/Ef remains approximately constant. At p < p n< j, 
d\n Ef /dlnp — 1/4. But for p > p n< j, it becomes smaller, 
and within the capture layer itself, extremely small, and 
the capture layer becomes correspondingly thick (see bot- 
tom panel of Figure 5). In fact, around neutron drip, the in- 
creased width of capture layers, coupled with the increasing 
number of capture layers per unit depth, makes the capture 
layers overlap. 

Because we only integrate one capture reaction at a 
time, our code cannot deal with overlapping layers, and 
so we artificially disregard several reactions indicated in 
Haensel & Zdunik (1990b). This should not lead to any se- 
rious problems, as our equation of state at these depths is 
insensitive to the exact (A,Z). Second, as shown in § 4.1, 
the quadrupole moments due to different capture layers add 
linearly, so overlapping capture layers can in principle be 
dealt with using superposition. 

2.3 Hydrostatic Structure of the Accreted Crust 

With the composition set, we now discuss the equa- 
tion of state. The electron pressure is that of fully de- 
generate relativistic T — particles, p e = 5.77 x 
10 29 erg cm -3 ( Pll /p c ) 4/:i , where pn = p/10 11 gem -3 . Free 
neutron pressure is given by the p„(n„) fit of Negele & 
Vauthcrin (1973), where n n = X-np/m^. We neglect the 
ion pressure and the non-ideal and thermal electron ef- 
fects on the equation of state. The total pressure is then 
p(p) — Pe(p) + Pn(p)- The resulting p — p relation is shown 
in Figure 6, where we have plotted it in such a way as to 
exhibit the changing balance between electron pressure and 
neutron pressure. It agrees quite well with the more sophis- 
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Figure 6. Equation of state used in this work. For densities below 
neutron drip, the EOS is just that of relativistic electrons. Neu- 
tron drip is indicated by a dashed vertical line. We have plotted 
the EOS in a manner such that the y-axis is p e ~ 4 / 3 for p < p n( j. 
The large steps at p < p n( j are from the electron captures, whereas 
at densities above neutron drip, the density contrasts weaken. 



ticated treatments of Haensel & Zdunik (1990a) and Brown 
(2000). 

In order to construct the crust, we pick a starting radius, 
pressure, and mass, and integrate the Newtonian equations 
of hydrostatic balance and mass conservation, 



dp 
dr 



-pg, 



dM r 
dr 



47rr 2 p, 



(18) 



(where M r is the mass enclosed inside r, and g = GM r /r 2 ) 
together with the rate equation (17) for the appropriate 
species, down towards the core of the star.^ Following Brown 
(2000), we decouple the calculation of the thermal structure 
from the hydrostatic calculation. This is justified since the 
equation of state is nearly independent of temperature. Only 
the electron capture rate is temperature-sensitive, and so the 
absolute locations of the capture layers that we find could 
be in error by a few meters. However, since the effect we are 
studying depends only on relative motion of the layers, this 
inaccuracy is of no concern. We start the integration well 
above the crust, at p ~ 10 21 erg cm -3 , in order to later ap- 
ply the thermal boundary condition. The crust begins where 
the ratio of Coulomb energy to thermal energy, 



t We initially used a Runge-Kutta integrator with adaptive step- 
size control. However, the stepsize adjustment algorithm tended 
to take steps that proved too large when we used them to solve the 
thermal structure and elastic perturbation equations. Namely, the 
very uneven mesh generated by the integrator led to large round- 
off errors and convergence difficulties. Our practical solution was 
to limit the maximum step size to a small fraction (10 — 6 ) of the 
radius. This generates a mesh that is mostly even and has extra 
resolution near capture layers as necessary. 



(where a — (3/47rn) 1//3 is the internuclear spacing) exceeds 
170. We stop the integration when we have reached the fidu- 
cial density p « 2 x 10 14 gcm -3 at the crust-core boundary 
(for a review and recent results see Pethick et al. 1995). Our 
(purely Newtonian) fiducial crust is 1.1 km thick and has a 
mass of 0.06Mq, with a mass of 1O _4 M0 in the outer crust. 
Brown (2000) points out that at high enough accretion rates, 
the combination of the low nuclear charge Z characteristic 
of the accreted crust and high temperatures may melt the 
crust (i.e., make Fcoui < 170) around neutron drip. We take 
this possibility into account by rerunning our Q22 calcula- 
tions with the top of the crust at neutron drip. 



2.4 Steady-State Thermal Structure of the Crust 

With the hydrostatic structure in hand, we compute the 
steady state thermal profile of the crust by solving the heat 
equation (without GR correction terms) 



V • F = pe = p (e nuc - e„) , 
where the thermal flux F obeys 
F = -KVT . 



(20) 



(21) 



Here e nuc is the local energy deposited by nuclear reactions, 
e v is the local energy loss due to neutrino emission, and K 
is the thermal conductivity. We neglect the compressional 
heating, which is negligible compared to the nuclear energy 
release (Brown 2000). In our calculation we mostly follow 
Brown (2000) and Brown &: Bildsten (1998) as far as the mi- 
crophysics is concerned, except for our treatment of nuclear 
energy release due to reactions in the crust. In particular, 
for e v we adopt a formula based on liquid phase electron 
vv bremsstrahlung (Haensel et al. 1996, Eq. [8]) because we 
expect the crust to be quite impure and hence neutrino emis- 
sion due to electron-ion and electron-impurity scattering will 
dominate over phonon scattering. For the same reason we 
use the electron-ion scattering conductivity (Schatz et al. 
(1999), Appendix A, generalization of Yakovlev & Urpin 
(1980) results) when computing K, rather than phonon- 
mediated scattering. 

Brown (2000) and Brown & Bildsten (1998) approxi- 
mated the heat deposition due to nuclear reactions as being 
uniform in the region around neutron drip. In this paper, 
since we are resolving individual capture layers, we also re- 
solve the heat release from them (though we find that this 
more accurate treatment does not lead to significant differ- 
ences). Let E nuc be the energy (ergs/nucleon) deposited in 
the transition layer by a single reaction (usually a pair of 
e - captures, accompanied by neutron emission at p > p n d)- 
The energy generation rate in the transition layer is then 



P^nuc — 



1 dX 



X t dr \4ivr 



(- 



m b 



(22) 



where X t is the mass fraction of the source nucleus at the 
top of the layer. The total energy (erg s -1 ) deposited in the 
transition layer is E nuc M /rrih- 

Equation (20) requires two boundary conditions. The 
boundary condition at the bottom of the crust is set by the 
ability of the core to radiate the heat flux from the crust as 
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Figure 7. Left (a): Thermal structure of Model S (0.5MEdd accretion rate, superfluid gap A = 1 MeV). Top panel shows the temperature 
of the crust. Middle panel shows the heat flux F r in the crust (F r < means the heat is flowing toward the core). Bottom panel shows 
the local energy release pe nuc (in erg cm" 3 due to reactions in the crust (sharp spikes) and the local rate of neutrino cooling pe v 

(also in erg cm" 3 s" 1 ). Right (b): Same as (a), but for Model N (0.5MEdd accretion rate, superfluid gap A = MeV, i.e. normal fluid 
in the core). 



neutrinos. For the bottom boundary condition we consider 
two cases: a normal core and a superfluid core. For the nor- 
mal core we use modified Urea neutrino emissivity (Shapiro 
& Teukolsky 1983, equation [11.5.24]). If the core is in a 
superfluid state, neutrino emissivity is suppressed^ by a fac- 
tor exp(— A/fceT), where A is the superfluid gap energy. In 
general, the gap energy will vary with the density in the 
core. However, we approximate superfluid effects by just in- 
cluding an overall exponentional factor in the modified Urea 
formula, 



5.3 xlO 39 erg s _1 



M 



Pnuc 

p 



1/3 



T98eXP (^)' (23) 



t We neglect neutrino emission via e-e bremsstrahlung, which is 
not exponentially suppressed by nucleon superfluidity and dom- 
inates neutrino emission from superfluid NS cores at T < 10 8 K 
(Kaminker & Haensel 1999) 



where p nuc is the nuclear density. If the core is composed 
solely of normal particles then A = 0, while the typical 
values of gap energy for superfluid cores is A w lMeV. The 
boundary condition at the bottom of the crust is just that 
all heat going into the core comes out as neutrinos, F + 
L COTC /4-Kr 2 = 0, where F is the radial heat flux; F = Ff. 

When nuclear burning in the upper atmosphere is 
steady, the outer boundary condition is set by the tempera- 
ture at the hydrogen/ helium burning layer, roughly 



T burn ^5.3x 10 8 k(^—) 2/7 , 

\m,Edd/ 



(24) 



where rn is the local accretion rate (Schatz et al. 1999). 
At sub-Eddington accretion rates (rh < (0.1 — l)rhEdd, the 
exact boundary is not known) the burning in the upper at- 
mosphere is not stable, leading to type I X-ray bursts (see 
Bildsten 1998b for a recent review). In that case the outer 
boundary condition is more complicated. However, Brown 
(2000) found that for high accretion rates, the local heating 
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in the deep crust makes the temperature there very insensi- 
tive to the outer boundary temperature. That is no longer 
true for low accretion rates, where the temperature in the 
inner crust is much more sensitive to the outer boundary 
temperature. While the thermal time at the burning layer 
is quite short, the thermal time in the majority of the crust 
is on the order of years, much longer than either the burst 
duration (~ 10 s) or the burst recurrence time (hours to 
days). Hence, we expect that in the time-averaged sense, 
the outer boundary temperature approaches that given by 
the steady calculation. Thus we adopt a very simple bound- 
ary condition by setting T = Tbu rn at p = 10 21 erg cm~ 3 
(the approximate location of the hydrogen/helium burning 
layer). 

Because superfluidity of the core changes the flux pro- 
file in the crust, we ran our calculations for two models. 
Both have the same accretion rate M = 0.5M B dd, use identi- 
cal microphysics, and share the same hydrostatic structure. 
However, model S (Figure 7a) has a superfluid core, with 
gap energy A = 1 MeV, while model N (Figure 7b) has a 
normal core (i.e. A = 0). These models are very similar to 
the ones obtained by Brown (2000), and we refer the reader 
to that paper for an in-depth review and discussion of their 
overall thermal properties. 



3 TEMPERATURE PERTURBATIONS DUE 
TO PHYSICAL ASYMMETRIES 

As discussed in § 1.3, there are several possible causes of 
spin-misaligned lateral temperature variations in the crust. 
Rather than simply assume a given temperature contrast 
ST (from which we can calculate Q22), in this section we 
calculate just how large either the composition variations 
or nuclear heating variations must be to imprint a given 
ST. This tells us how large either of these effects must be 
to generate a sufficient Q22 for gravitational radiation to 
balance the accretion torque. 

3.1 Possible Sources for the Temperature 
Variations 

Lateral differences in the crustal composition will have two 
effects. The first is lateral variations in the amount of en- 
ergy deposited in the crust by nuclear reactions (since nu- 
clear transmutations of different elements deposit different 
amounts of energy) . We denote lateral variations of this type 
as / nuc = SE nuc /E DUC , so that a nonzero /nuc means that 
more energy is released on one side of the crust than the 
other§. Second, the charge-to-mass ratio Z 2 /A is likely to 
differ if burning proceeds to different A on different sides 
of the star. The conductivity in the crust scales as K oc 
(Z 2 /A)~ 1 pT nk (Schatz et al. 1999), while neutrino emissiv- 
ity is e„ oc (Z 2 /A)pT n ° (Haensel et al. 1996). For the micro- 
physics employed here, n k ~ 1 and n e « 6, so that the con- 

§ Since the rate at which the energy released in the crustal nu- 
clear reactions depends on the local accretion rate (see Eq. [22]), 
a similar effect would occur if the crust had laterally uniform 
composition entering at the top, but different local compression 
rates. 



ductivity and neutrino emissivity will vary laterally with the 
composition variation. We denote the Eulerian perturbation 
in the charge-to-mass ratio as / CO mp = 8{Z 2 / A) / {Z 2 / A). Re- 
gions with /comp > are more opaque and radiate neutrinos 
more efficiently. The hydrostatic structure is hardly affected 
however, as the dependence of the EOS on A and Z is very 
weak. 

Both / nuc and / mmp will lead to a calculable lateral 
temperature variation, ST. An important issue to clarify is 
the role of the core, which has high thermal conductivity 
and hence is nearly isothermal. Although the core tends to 
smooth out temperature variations (i.e., decrease ST), we 
show that it does not force them to zero in the crust. Be- 
cause the core's thermal conductivity is much higher than 
the crust's, ST in the core is much smaller than in the crust. 
Hence, we approximate the core as perfectly conducting and 
isothermal. The crust is internally heated by nuclear reac- 
tions near neutron drip and therefore is not at the same 
temperature as the core. The crustal thermal equilibrium is 
set mostly by the heat flux in the radial direction, implying 
that the radial temperature gradient on the side with posi- 
tive /nuc or /comp must be steeper in order to connect to the 
same core temperature. We now calculate this ST. 

3.2 The Thermal Perturbation Equations and 
Boundary Conditions 

While the ultimate effect of / nuc and / CO mp is to shift the 
location of capture layers, to first order we can calculate ST 
by considering the effects of / nuc and / mmp on a spherically 
symmetric background. After computing ST, we will find in 
§ 4 how the crust hydrostatically readjusts. To first order, we 
then set 8p = 0, and the conductivity perturbation equation 
is simply 

SK _ ST 

— ~ /comp + flfc-y", (25) 

while the neutrino emissivity perturbation is = 
/comp + n e 8T/T. 

The nuclear energy generation rate perturbation re- 
quires special attention. Nuclear reactions are generally 
quite temperature-sensitive. However, despite the temper- 
ature perturbation, the total energy release of a complete 
capture layer depends only on the local accretion rate and 
the total E nuc of the element. A temperature perturbation 
shifts the capture layer and hence leads to a local change 
in e nU c on scales smaller than the distance over which the 
layers shift. However, it cannot change the total amount of 
energy released. Since we neglect the shifts of capture lay- 
ers at this stage, we simply have 5e n uc/enuc = /nuc- In some 
sense, we average the energy generation rate over the scale 
of the entire capture layer. This approximation simplifies 
the calculation considerably. The Eulerian perturbation to 
e = e nU c — £k (the local nuclear heating rate minus the neu- 
trino cooling rate) is therefore 

8e _ gnuc ( f 8T\e v 

— — /nuc ~ I Jcomp "T I — • \^®) 

We assume that the angular dependence of all per- 
turbed quantities is oc Yi m , i.e. 8T(r, 6, 4>) — ST(r)Yi m (8, <f>). 
Perturbing the heat equation (20) and the flux equation (21) 
and keeping only first-order terms gives 
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(SK dST \ ^<"Vv 



(27) 



and 

re - "(tt-t) 

+ 4(f)- <-> 

Substituting for 8K/K and 5e/e, and using Eq. (21), 
Eq. (28) becomes 



1 d / 2R ,rf5T\ F ^dST in ,..Ki 



K5T 



F 2 I ST 

- \ n kJ^f +p(n fe (e nU c -e„) + n e e„) ^ — 



Jcomp (2£;y £nuc) /nuc^nuc 



(29) 



where we have neglected drik/dr, and for simplicity have 
taken the composition perturbation to be radially uniform, 
so dfcomp/dr = 0. 

The thermal perturbation problem, Eq. (29), requires 
two boundary conditions. At the top of the crust, the exact 
boundary condition can be obtained by matching to a flux- 
temperature relation in the ocean, where the thermal profile 
is set by compressional heating and that portion of the nu- 
clear energy released in the deep crust that flows upward 
through the ocean, rather than down into the core (Bild- 
sten & Cutler 1995; Brown & Bildsten 1998). But since 
the ocean's thermal conductivity is much higher than the 
crust's, we simplify our calculation by adopting the bound- 
ary condition, 5T(top) = 0. The variation in the flux coming 
out of the crust is not zero; in fact it is potentially observ- 
able. 

Now consider the boundary condition at the crust-core 
interface. The thermal conductivity in the core is at least 
several orders of magnitude higher than in the crust. Hence, 
any extra flux into the core can be carried with a very small 
temperature perturbation, as we now show. First consider a 
superfluid core. The boundary condition for the spherically 
symmetric calculation is F CO r C = 0, since neutrino emission is 
suppressed and the core cannot radiate away any significant 
heat flowing into it (§ 2). However, we are now lifting the 
restriction of spherical symmetry, so heat can flow into the 
core on one side and out the other. The size of ST in the core 
is then related to the magnitude of the radial flux perturba- 
tion by <5F r , corc ~ (K COIC T/R)(ST/T)\ COIC . The magnitude 
of the transverse flux 5F± tCOie is of the same order, contrary 
to the situation in the crust, where the transverse heat flux 
is much smaller. The radial flux perturbation in the crust 
is 5F r , crust ~ {K CIUSt T/AR)(ST/T)\ CIUSt , where (ST/T)\ CIUSt 
is the typical magnitude of the temperature perturbation in 
the crust, and AR is the thickness of the crust. Continuity 
of the radial flux at the crust-core boundary then gives 



ST 
T 



R 8T_ 
"AR ~T 



< 



ST 
~T 



(30) 



i K C OYC 

When the NS core is not a superfluid, it can emit neutri- 
nos and the equilibrium model has a nonzero flux F coro 
going into the core (see Figure 7). The radial flux per- 
turbation in the core is then determined by the com- 
petition of two terms in Eq. (27), (SK/K)F COTC and 
KdST/dr ~ (K coie T/R)(ST/T). The core is nearly isother- 



mal, so K COIC T/R 3> -Fcoro (i.e., the proper estimate of 
Fcoro is K COIC AT/R, where AT <§C T is the difference in 
temperature between, say, the center and the crust-core 
boundary). Therefore, the second term, KdST/dr, domi- 
nates, and, just as in the case of a superfluid core, we have 
SF TtCOTe ~ (A' CO r c r/_R)((5r/r)|coro. Thus we again arrive at 
Eq. (30), i.e., the typical magnitude of ST in the crust is 
several orders of magnitude larger than ST in the core. 

Hence for both the normal and superfluid core, to good 
accuracy we can take ST = as our boundary condition 
at the crust-core boundary. With this boundary condition, 
we do not need to model the core, as we effectively assume 
that it is perfectly conducting. However, perturbed flux is 
flowing through the core. 



3.3 The Resulting Temperature Variations 

Solutions of the perturbation problem (29) are shown in 
Figure 8 for the quadrupole (I = m = 2) case^, for ac- 
cretion rate M = 0.5MEdd- In Figure 8a we presume that 
the nuclear heating varies laterally by 10%, i.e. / nuc = 0.1, 
but take / CO mp = 0, so conductivity and neutrino emissivity 
are unperturbed. The resulting temperature perturbations 
are shown in the top panel, with the solid line represent- 
ing the superfluid-core case (gap energy A = 1 MeV) , while 
the dashed line corresponds to the case of a normal core 
(A = 0)JI Figure 8b shows the solutions for the opposite 
case, where / nuc = but / comp = 0.1. Bottom panels of the 
figures show the perturbation in the radial heat flux SF r . 
The transverse flux in the crust is negligible in all cases, and 
is not shown here. 

We see that the typical ST /T for M = 0.5M Edd near 
the bottom of the crust is roughly (0.1 — 0.3)/, where / is 
either / CO m P or / nu c- Equivalently, the temperature pertur- 
bation ST is a few xlO 6 K for / mmp or / nuc of 10%, with 
the maximum ST attained around neutron drip. The tem- 
perature perturbations tend to be larger in models with a 
normal core than with a superfluid core. One might have 
expected that the amplitude of ST /T would be of order / nuc 
or /comp- This is not the case because the crust can radiate 
away some of the "extra" flux it needs to carry by emitting 
neutrinos, thus reducing the lateral temperature gradient. 
However, this enhanced neutrino emission does not com- 
pletely eliminate the lateral temperature gradient. 

In Figure 9, we survey the dependence of the tem- 
perature perturbations on the accretion rate. Most of the 
heat released by the nuclear reactions near neutron drip 
(roughly 1 MeV/accreted baryon) flows towards the core 
and is lost as neutrino emission either from the crust or 
from the core. However, a fraction of the heat flows to- 
wards the surface of the star. This flux F(r top ) is plotted 



^ The thermal perturbation equations do not depend on m, so 
our solutions are valid for any 1 = 2 perturbation. 
II We have likely overestimated the neutrino emission suppression 
in a superfluid core (i.e. our model S). As Brown (2000) points out, 
the neutrino emissivity depends very sensitively on temperature 
so that even a moderate number of normal particles in the core 
can radiate away an appreciable fraction of the flux from nuclear 
reactions in the deep crust. Hence it is likely that ST/T of model N 
is more representative. 
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Figure 8. Left (a): Temperature and flux perturbations in the crust resulting from laterally varying energy release, / nuc = 0.1, but 
unperturbed conductivity, / comp = 0. Angular dependence of the perturbation is I = 2, i.e., a quadrupole. Top panel shows the 
temperature perturbation 8T/T, while the bottom panel displays the radial flux perturbation SF r . Solid lines correspond to Model S 
(superfluid core, M = 0.5MEdd)> while dashed lines represent Model N (normal core, same accretion rate). Transverse heat flux 8F± is 
not shown, as it is much smaller than SF r in the crust. Right (b): same as (a), but the energy release is not perturbed, / nuc = 0, while 
the charge to mass ratio Z 2 /A is varied laterally, / mmp = 0.1. 



in the top panel of Figure 9, scaled by the accretion flux 
F A « (1/4tyR 2 )(GMM/R), or w 200 MeV/47ri? 2 /accreted 
baryon. The solid line denotes a model with a normal core, 
while the dashed line corresponds to the model with a su- 
perfluid core. In the superfluid case a smaller fraction of the 
heat escapes to the surface because the crust is hotter and 
hence crustal bremsstrahlung is more effective at radiating 
the heat deposited by nuclear reactions on the spot. 

The middle panel of Figure 9 shows the magnitude of 
the temperature perturbations that result when either the 
local heating rate is perturbed by 10% (/ nu c = 0.1, lines 
marked with filled triangles) or the charge-to-mass ratio is 
altered, leading to lateral variations in the conductivity and 
neutrino emissivity (/ C om P = 0.1, lines marked with open 
circles). The typical magnitude of 8T/T is a few percent for 
the models with a normal core (dashed lines) and is some- 
what smaller for models with a superfluid core (solid lines) . 
We must emphasize that our temperature perturbation cal- 
culation is linear, and hence a 20% asymmetry will result in 
temperature perturbations that are a factor of 2 larger, so 
long as ST/T < 1. 

These temperature perturbations displace the capture 
layers and create a quadrupole moment. But in addition, 
they result in lateral variations in the flux leaving the top of 
the crust. The amplitude of this flux perturbation, SF r (r t0 p) 



is plotted in the bottom panel of Figure 9, scaled by the 
accretion flux Fa- The typical magnitude of the flux modu- 
lation is 8F r /FA ~ 10~ 4 (//0.1). These hot and cold spots, 
when moving in and out of the view of the observer due to 
the rotation of the neutron star, can generate a modulation 
in the persistent emission. Observational implications of this 
effect are discussed in § 5. 

The existence of lateral temperature variations depends 
on the continual heating of the crust. When accretion halts, 
thermal diffusion with slowly equalize the temperature lat- 
erally. How long will this process take? 

First, let us compare the transverse flux perturbation 
in the crust SF± to the radial one, 5F r , as one would ex- 
pect that lateral heat transport could be responsible for 
equalizing the temperature. From (21), the background flux 
in the crust is F ~ KT/AR, while from (27) the trans- 
verse flux perturbation SF ± = KST/R ~ F(AR/R)(5T/T). 
The perturbed radial flux is the sum of two terms, both of 
size ~ F(ST/T). So for low angular order I, the perturbed 
transverse flux is smaller than the radial piece by a factor 
~ AR/R. Hence, the transverse heat flux in the crust does 
not wash out the temperature perturbations. The heat is 
transported predominantly radially through the crust, and 
both radially and transversely through the core. Hence, if ac- 
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Figure 9. Top panel: The unperturbed flux emerging from the 
top of the crust due to the (spherically symmetric) direct heat 
input by the nuclear reactions around neutron drip. The flux 
F(rtop) has been scaled by the accretion flux F\ = 1.4 X 
10 25 erg cm 2 s _1 (M/M Edd ), so F/F A = 10~ 3 corresponds to 
a luminosity of 200 keV per accreted baryon emerging from the 
top of the crust. Solid line denotes the model with a supcrfluid 
core (A = 1 MeV), dashed line denotes the model with a nor- 
mal core (A = 0). Middle panel: the maximum magnitude of the 
crustal temperature perturbation ST/T induced by asymmetric 
heat sources (/nuc = 0.1, lines marked with filled triangles) and 
opacity/neutrino emissivity variations (/ CO mp = 0.1, lines marked 
with open circles). Solid lines are the results for the model with 
a superfluid core, dashed lines correspond to the model with a 
normal core. These results scale linearly with / nuc and / CO mp- 
Bottom panel: the perturbation to the flux emerging from the 
top of the crust due to the asymmetry in nuclear energy release 
(/ nuc = 0.1) and Z 2 /A (/ mmp = 0.1). The flux perturbation SF r 
has been scaled by the accretion flux Fa. The legend is same as 
in the middle panel. 



cretions stops, the temperature perturbation will be washed 
out in a thermal time at the bottom of the crust, 



*th ~ 6 yr 



P 



10 32 erg cm -3 



3/4 



(31) 



(Brown et al. 1998). The transverse temperature asymmetry 
will therefore persist so long as (1) accretion persists on 



timescales longer than a few years, and (2) the crust has 
compositional asymmetries. 

What does this mean for our quadrupole generation 
mechanism? In steadily accreting sources, we would expect 
the quadrupole moment not to vary in time, and hence their 
spin will be set by the competition of accretion with grav- 
itational radiation torque. However, as suggested by Bild- 
sten (1998a), transiently accreting sources with long recur- 
rence times may be able to spin up to shorter periods, since, 
because of the short thermal time, the temperature varia- 
tions in the deep crust and their quadrupole moments may 
be lower than what one would expect from just the time- 
averaged accretion rate. 



4 THE ELASTIC DEFORMATION OF THE 
CRUST 

In this section we derive and solve the perturbation equa- 
tions that describe the elastic response of the crust to lat- 
eral composition gradients. Though our primary interest is 
rotating NSs with v a ~ 300 Hz, we consider deformations 
of nonrotating, spherically symmetric background models. 
For static deformations (i.e., no Coriolis force), centrifugal 
terms modify our results only by order (v 3 /v b ) 2 ~ 5%, where 
~ 1-5 kHz is the NS breakup frequency. A given level of 
temperature or composition asymmetry gives rise to "moun- 
tains" of a certain size, where that size is only slightly modi- 
fied by rotation. However, because the deformed star rotates, 
it emits gravitational waves. 

We treat two cases: one where the lateral composition 
gradient is due to a wavy e~ capture boundary (and the 
gradient is therefore confined to a region near the layer), and 
another where the lateral composition gradient is uniform 
throughout the crust. Real neutron stars probably exhibit 
both types of composition gradients at some level; for small 
distortions, their effects should add linearly. These lateral 
asymmetries, treated as perturbations on a homogeneous 
background, cause pressure imbalances that source a crustal 
displacement field £ a . Treating the crust as an elastic solid, 
we solve for the £ a that brings the crust back to equilibrium, 
with the gravitational, pressure, and shear-stress forces all 
in balance. We then compute the density perturbation 8p 
and the resulting Q22- An important underlying assumption 
is that the stresses are small enough that the NS crust can 
be in static equilibrium. More precisely, we assume that the 
crust is deforming slowly enough that the elastic part of the 
viscoelastic stress tensor dominates over the viscous part. 
If the stresses are much larger than the crust's yield stress, 
then this cannot be true. 

The process by which deformation is built up in a NS 
is undoubtedly complex, involving the viscoelastic response 
of the crust to temperature and composition gradients that 
are built up over time, during which the primordial crust is 
replaced by accreted matter. A full understanding of crustal 
deformation would involve solving the time-evolution equa- 
tions for all relevant aspects of the crust, from the moment 
accretion starts. This would involve plastic flow or break- 
ing whenever the crustal yield strain was reached. Such a 
"movie" of the crust's history is well beyond our current 
abilities. Instead, our solution of the perturbed hydro-elastic 
equations essentially amounts to taking a "snapshot" of the 
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crust. We find that this approach gives a lot of detailed in- 
formation, but necessarily we must put in some things "by 
hand." 

One important thing we put in by hand is the reference 
state of the crust: for simplicity, we take it to be spherically 
symmetric. Namely, we assume that if one could somehow 
"undo" the "extra" e - captures that cause the capture layer 
to be wavy, then the neutron star would "bounce back" to a 
configuration with zero quadrupole moment. Since the neu- 
tron star crust may have undergone a long history of plastic 
flow and cracking, it is not clear how closely the real crust 
matches this picture. Fortunately, it is easy to see how this 
assumption affects our results. Assume that, in the absence 
of any lateral composition gradients currently driving the 
neutron star away from spherical symmetry, the star would 
adopt a shape with multipole moments Qi^- The super- 
script "hist" refers to the fact that for "historical reasons," 
the crust has evolved to an equilibrium shape that is non- 
spherical. (Of course, almost all the solid objects we use 
in daily life have nonspherical equilibrium shape: we man- 
ufacture them that way.) Let Qf^ be the piece due to the 
current lateral composition gradient, assuming the reference 
shape is spherical. Both deviations will be small, and to first 
order they add linearly: 



Q total 
illl — 



■jpcrt 
' Im 



(32) 



We guess that the dominant historical forces that shaped 
Q22 St ( e -g-, through viscoelastic flow) are the same ones that 
currently shape Q£| r . Correspondingly, we guess that the 
two pieces tend to add coherently, rather than to cancel. 
For example, consider a spherical shell of steel, and apply a 
large inward force at two opposite points on the equator. At 
first, Qf^ al is j us t Qfm ' tne distortion due to the existing 
force. But over time the steel also relaxes somewhat, and 
the effect is obviously to increase the total deformation. One 
might then wonder how large Q^ 3,1 could grow over time, 
for a fixed Q^*- ^ i s partly to address that question that 
we derive, in § 6, an upper limit limit on Q^ 3 " 1 , independent 
of the relative contributions of Q^* ana Qfm'- 

In the rest of § 4 we omit the superscript "pert" from 
Q2T ) but hopefully it is now clear that the neutron star's 
total Q22 does have another piece, Q 2 2 St - 

Now let us turn to the source of Q^"* '■ As stated above, 
we consider two types: composition gradients due to wavy 
capture boundaries, and composition gradients that are ra- 
dially uniform. Our model of the wavy capture boundary is 
as follows. We posit the existence of some Eulerian tempera- 
ture perturbation ST = Re{8T(r)Yi m (8, 4>)}, superimposed 
on a background that would otherwise have spherically sym- 
metric composition. One possible source of such ST is the 
asymmetric heat flow due to laterally varying nuclear heat- 
ing or conductivity (Z 2 /A), as computed in § 3. Regions 
where 5T(r) is positive (negative) have their e~ capture 
layers shifted to lower (higher) density, as discussed in § 1.2 
and § 2.1. Essentially one ends up with a crustal EOS that is 
angle-dependent near the capture layer, and which requires 
shear stresses to be in equilibrium. 

The other source we consider is a (radially) uniform 
lateral composition gradient, Ap c / p c = Re{C Yi m (9, </>)}, 
for some constant C. The resulting Q22 will scale linearly 
with C; our fiducial choice is C — 5 x 1CF 3 . The effect is 
again to give angular dependence to the crust's EOS, but 



this dependence is now small and uniform rather than large 
and confined to the e~ capture regions. It should be clear 
that we do not think Ap e /p e (r) will really be a constant, 
but we just consider it in an averaged sense. Our use of 
Lagrangian Ap e here instead of Eulerian 5p c reflects the 
fact that we are specifying the composition gradient that 
would exist (due to accretion and/or burning asymmetries, 
say) if the crust did not elastically adjust. 

Having explained the background solution and the 
sources, we now derive the crustal perturbation equations. 



4.1 The Crustal Perturbation Equations 

Our calculation is based on the following assumptions and 
approximations. We use Newtonian gravity throughout. We 
model the crust as an elastic solid with shear modulus p and 
two-parameter equation of state: p = p(p, p e ), where p is the 
density and p c is the electron mean molecular weight. We 
neglect temperature in the equation of state, except insofar 
as it affects p a (see Appendix A for details). We neglect 
the slight overall downward flow of matter due to accretion; 
this is justified since the ram pressure at relevant depths is 
completely negligible compared to gravity or shear forces. 
We presume that the crust responds elastically, i.e., we do 
not allow the stresses to relax plastically or via cracking. 
The stress-energy tensor of the solid is then 



Tab 



-P(P, Me)ffa6 + H [ V £& + V 6 £ a - ^Se^V^c ) , (33) 



where g a is the flat 3-metric and V a is its associated deriva- 
tive operator. The equation of static balance is V°t & = 
pV 3> where <E> is the gravitational potential. 

We carry out first-order perturbation theory with a 
spherically symmetric background model (constructed in 
§ 2) that has zero shear stress. Thus we treat £° as a first 
order quantity. We neglect the perturbation of the gravi- 
tational potential $ (the Cowling approximation; but see 
§ 6.5 where we relax this approximation). We use S to rep- 
resent Eulerian perturbations and A for Lagrangian per- 
turbations. For scalar quantities A, they are related by 
AA = 8A + £ a V a A, where £ a is the displacement vector of 
the elastic solid from its original state. Since the background 
model is spherically symmetric, both the Eulerian and La- 
grangian variations of scalars are proportional to Yj m , and 
AA(r) = SA(r) + £ r dA/dr. 

The Lagrangian pressure perturbation depends on 
which effect causes the crustal EOS to have non-trivial an- 
gular dependence. In the case of a smooth composition gra- 
dient, Ap e /p e , we write 



Ap 



dp 
dp~ 



Ap + 



dp 
dp c 



Ap e 



(34) 



since p = p(p, p c ). For the ST source term, we write 

AT (35) 
dT^ 



A ^ 

dp 



a dp 
A P+ JT- 



dp 



Ap + 



dp c 



where in the second equality we used the relation between 
Eulerian and Lagrangian perturbations. Evaluation of the 
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source terms and coefficients in equations (34) and (35) is 
discussed in detaii in Appendix A. The oniy transformation 
we make here is to decompose the source terms into spherical 
harmonics, i.e., we write A/x e = Ap e (r)Yi m (6 ,<f>) and ST — 
5T{r)Yi m (9,(j)). In either case, the perturbation "sources" a 
displacement vector £ a of the form 



C = Z r (r)Y lm r a + ^l3-\V a Yi m , 



(36) 



where f a is the unit radial vector and (3 = \Jl(l + 1). 
The perturbation equations for the crust are then** 

V a 8r ab = 8pgf b . (37) 

Here g(r) = GM r /r 2 is the local gravitational acceleration 
of the background model, and 8r ab is given bytt 

Sr ab = gabYimSTrr + e a bYi m [2n(£ T /r — d T £r)] (38) 
+ fabSr r ± + A ab 2p/3^±/r, 

where 



Sr r ± 
and 



-6p + lt[-drZT ~ jfr/r+^l/f] 

rdr(£,±/r)+fc ± /r 



e ab = g ab - r a r b , 

fab = P~ l r{f a VbYi m +rbV a Y lm ), 

A ab = r 2 r 2 V a V b Y lm + f ab . 



(39a) 
(39b) 



(40a) 
(40b) 
(40c) 



Note that there is no Sp, term in Eq. (37), since in Eq. (33) 
the shear modulus p is multiplied by terms involving £ a , 
which vanish in the background model. The Eulerian den- 
sity perturbation Sp in Eq. (37) follows from the perturbed 
continuity equation, 

S P = -V a ( P ia) = -{±£ (r 2 p&) ~ P^} Y lm , (41) 

and the Eulerian pressure perturbation Sp — Ap + pg£ r in 
Eq. (39a) follows from either Eq. (34) or Eq. (35), depending 
on the type of perturbation imposed upon the crust. 

The radial and perpendicular pieces of Eq. (37) yield 
two second-order equations for the variables £ r and £±. Fol- 
lowing McDermott et al. (1988), we rewrite these as four 
first-order equations, using variables 



21 



23 = 



k 

r ' 

£± 
f3r' 



22 



2 4 = 
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Ar r± 


8T r± 


Pp 


Pp ' 
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dlnp 
dlnr 



(42) 



dzi 
dlnr 



1 + 2 21 H 2 2 

V «3 a 3 J as 



(43a) 



** These equations can also be derived by explicitly writing out 

the components of the stress tensor and its divergence in spherical 

coordinates. We found the tabulation of the components of the 

stress tensor in Takcuchi & Sato (1972) very useful. 

tt When not explicitly specified, all physical quantities are 

obtained by taking real parts of the corresponding complex 

expressions. 



+ 1(1 + 1)— 23 + — AS, 

Q3 Q3 



d In r 



dZ2 = (UV~W + 12T^-4^±) Z1 (43b) 

V Q3 Q3 / 



23 



+ J(J + 1)Z4 - 4^AS, 

a 3 

dz z 1 

-T, = -21 H 2 4 , 

amr ai 
dz4 
dlnr 



(43c) 
(43d) 



= I V — 6r h 2 21 22 

V «3 «3 / «3 

+ — { [21(1 + 1) - 1] aia 2 + 2 [1(1 + 1) - 1] aj) z z 

Q3 



+ (V-3)2 4 + 2^AS, 



where 



d In o 
(7^—^+2, 
dlnr 



V = 



pgr 

P 



dlnp 
dlnr ' 



(44) 



/ „ 2ii „ 4« 

ai=n/p, a 2 =r---, a 3 = r+--. 

The terms T, 04, and the source term A5 depend on the 
type of perturbation. In the case where the perturbations are 
due to a lateral composition gradient Ap e , the derivatives 
in Eq. (34) are carried out at constant composition p c , and 
we have 



r = 



<91np 



<91np 



a 4 = 0, AS = 



dlnp 



din p c 



Ap c 

He 



(45) 



On the other hand, in case of the ST source term, the per- 
turbations in Eq. (35) are at constant temperature, so 



r = 



dlnp 



dlnp 



a 4 



d Inp 



dlnT 



dlnT AO dlnp 
dlnr' <91nT 



ST 
~T 



(46) 



We describe how to compute the various thermodynamic 
derivatives in the above equations in Appendix A. 

Except for the terms involving a.4 and AS, which we 
have highlighted by writing them in boldface, Eqs. (43) are 
the same as the zero-frequency limit of Eqs. (13a-d) in Mc- 
Dermott et al. (1988). However, in the case of adiabatic pul- 
sations considered by McDermott et al. (1988), one writes 
Ap — (dp/ dp)\ s Ap, where the partial derivative is taken at 
constant entropy. In this case, there are no source terms AS, 
no terms involving the temperature gradient dT/dr, and T is 
just the adiabatic index I\ = (dlnp/dlnp) s . The AS source 
term in Eq. (43) arises from the fact that our perturbations 
(Eqs. [34] or [35]) involve an explicit change in composition 
as well as a change in density. In addition, our V term is de- 
fined differently from McDermott et al. (1988), and depends 
on the type of perturbation. Consequently our terms a 2 and 
03, defined in Eq. (44), are also different from McDermott 
et al. (1988). 

Because Eqs. (43) are just linear equations with an in- 
homogeneous forcing term AS proportional to Ap e or ST, 
it is clear that the resulting quadrupole moment is linear in 
A/x e or ST, respectively. For the ST perturbations, the coef- 
ficient (<91np/<91nT)|p in AS vanishes except near capture 
boundaries, and so the total Q22 for the NS is just the sum 
of the (322's generated by each capture layer individually. In 
§ 4.3 and 5 we are therefore justified in looking at solutions 
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for a single capture layer, since multiple capture layers can 
be dealt with by superposition. 

4.2 Boundary Conditions and Solution Methods 

The solid crust lies between a liquid ocean and a liquid core. 
In the liquid, ai = and Q2 = 03 = T. By integrating 
Eqs. (43) across the crust-ocean boundary at r = r top , and 
the crust-core boundary at r = rbot, it is easy to see that 
zi, 22, and 24 (i.e., the radial displacement £ r , and tractions 
Ar rr and Ar r ±) must be continuous when going from the 
crust to a liquid. On the other hand, it follows from Eq. (43c) 
that 23 (the transverse displacement is allowed to have 
an arbitrary discontinuity at r top and r bo t- These connection 
conditions, however, require the knowledge of £ r , Ar rr , and 
Aryx in the liquid parts of the neutron star. We now show 
how to use physical considerations to remove this require- 
ment. 

For static solutions with I / 0, the Eulerian pressure 
perturbation Sp(r)Yi m must be zero everywhere in the liquid 
parts of the star, since a non-zero Sp would lead to a horizon- 
tal pressure gradient in the fluid Sp(r)V a Yi m , which would 
(in the absence of counterbalancing shear stresses or per- 
turbations of the gravitational field) cause the fluid to flow. 
Hence, at a liquid-crust boundary, we have Ar rr (crust) = 
— Ap(liquid) = pg(, r - Therefore, at r top and rbot, the Eule- 
rian STrr = 0. Using this, along with the constraint that 5r r ± 
must vanish in the fluid (inviscid fluid cannot support shear 
stresses) one arrives at the following boundary conditions 
(in agreement with McDermott et al. (1988)): 



Z2 — Vzi, 24 = 0. 



(47) 



Note that the above boundary conditions do not require the 
knowledge of displacements and tractions in the liquid parts 
of the star. 

Equation (47) assumes that there is no density disconti- 
nuity at the solid-liquid boundary. However, at the interface 
between the crust and the ocean there may be a discontinu- 
ity (p s — pi)/ p ~ 5 x 10~ 5 (where p s and pi are densities of 
solid and liquid at the interface) arising from the latent heat 
of melting of the crust. There may also be a density jump 
of a few percent associated with a phase transition between 
the inner crust and the core of the neutron star (Pethick & 
Ravenhall 1995; Pethick et al. 1995). When there is a den- 
sity discontinuity, the boundary condition at the edge of the 
solid is modified to 



22 



V 2i, 



24 = 0. 



(48) 



Of course, Eq. (47) is just a special case of Eq. (48), where 
p 3 = pi at the interface. 

Thus we have a system of 4 ODE's with 4 boundary con- 
ditions (2 at each boundary) , with non-zero source terms in 
the interior region. The system can be solved in any number 
of standard ways. For example, one can choose the values 
of 21 and 23 at the top of the crust, compute 22 and 24 
from the boundary conditions, then integrate to the bottom 
of the crust (using a Runge-Kutta type algorithm) and cal- 
culate the residuals of the two boundary conditions there. 
Then one can use a 2-dimensional Newton's method to find 
appropriate starting values z\ and 23. In practice, we found 
this method to be unstable. On the other hand, we find that 



a Henyey-type relaxation algorithm (as described in Press 
et al. 1992) converges very reliably. However, roundoff errors 
can cause trouble if the relaxation mesh is not fine enough 
and uniform, as already noted. 

Given the solutions of Eqs. (43), the perturbed multi- 
pole moments Qi m are computed in terms of the 2, variables 
by using Eq. (1), 



Qi 



rtop 



( Zl + + l)z 3 )r l+2 pdr - [r' +3 p 2l ] 



rtop 



(49) 



where rbot and r top are the bottom and top of the crust, 
respectively. If there is a density discontinuity at the crust- 
core or the crust-ocean interface, then p in the boundary 
term of (49) is replaced by pi. Since the calculation of the 
quadrupole moments involves subtracting two large numbers 
to get a small number, it is useful to have an independent 
accuracy check. An alternate expression for the quadrupole 
moment is derived by finding 5p from Eq. (37), integrating 
it using Eq. (1), using the boundary conditions and carrying 
out substantial algebra to show 

Qi m = - f ^ l\l(l + l)z 4 -2a 1 (2^+l(l + l)z,\ 



+ (l + 4-U)(z 2 - Vzi) \r l+2 pdr 



(50) 



Unlike Eq. (49), this expression for the multipole moment 
remains unchanged even when there are density jumps at 
solid-liquid interfaces. In practice, we refine our solutions to 
Eqs. (43) until Qi m computed from Eqs. (49) and (50) differ 
fractionally by less than 10~ 6 . 

Since the crust is geometrically thin, one may wonder 
if accurate results may be obtained by adopting a plane- 
parallel approximation, i.e., neglecting non-sphericity of the 
crust. In this case, the radial momentum balance equation 
is 



ik8r r ± H — —&T rr — Spg, 
dr 



(51) 



where the lateral (x) dependence of all perturbed quantities 
is e tkx and k ~ + 1) /R is the transverse wave number. 
Integrating the above equation from rbot to r top and using 
the boundary condition Sr rr = 0, we get a plane-parallel 
estimate for the quadrupole moment, 



Q22 



R 



J 



ikSr r 



-dr = l(l + l)R 



Z4prdr. 



(52) 



Compare the above expression to Eq. (50). In the plane- 
parallel case, only the Sr r ± component of the shear stress 
contributes to the quadrupole moment. As it turns out, 
the 5r rr component of the shear stress is much larger (see 
§ 4.3). In other words, the crust holds up the majority of the 
quadrupole moment by stretching like a spring with a weight 
attached to it, rather than bending like a crossbar. In the 
plane-parallel case, this vertical motion of the crust is fully 
taken into account. However, because of the artificial sym- 
metry of Eq. (51), the dominant stress term does not enter 
in the plane parallel expression for the quadrupole moment, 
Eq. (52). Quadrupole moments computed using the plane- 
parallel approximation (Eq. [52]) underestimate the values 
computed using full spherical symmetry (i.e., from Eqs. [49] 
or [50]) by large factors (~ 10). 
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Figure 10. Response of the crust to a temperature perturbation 
due to an / nuc =0.1 asymmetry of the nuclear heating rate in 
the crust. The background model has a normal core (gap energy 
A = 0) and the accretion rate is 0.5MEdd- Only one shallow 
capture layer (Q = 23 MeV) is activated. Top two panels show 
vertical and horizontal Lagrangian displacements, £ r and £x> m 
meters, the two middle panels show the stresses, Ar rr and Ar r j_, 
in 10 29 ergcm — 3 . The bottom panel shows the negative of the 
source term, — (d lnp/8T) p 5T (see Eq. [46] and Appendix A). 



4.3 The Nature of the Solutions 

We now describe the solutions to our crustal perturbation 
equations, both for the case where the capture layers are 
deformed due to local temperature variations ST, and for 
the case where the perturbations are sourced by a radially 
uniform composition perturbation Afj, e /fj, e . 

Figures 10, 11, and 12 show the response of the crust to 



Figure 11. Same as Figure 10, but for a displacement of a 
capture layer at intermediate depth, Q = 42.4 MeV. 



a temperature perturbation ST(r)Re{Y22} acting on shallow 
(threshold energy Q = 23 MeV), medium (Q = 42.4 MeV) 
and deep (Q = 95 MeV) capture layers, respectively. The 
background model has a normal core, and the accretion rate 
is 0.5M E dd- The background thermal structure is shown in 
Figure 7, and 8T(r) is plotted in Figure 8. The top two pan- 
els show the vertical (£ r ) and horizontal (£x) Lagrangian 
displacements of the fluid elements in the crust, as defined 
in Eq. (36). The next two panels show the Lagrangian com- 
ponents of the perturbed stress tensor, Ar rr and Ar r ±, de- 
fined in Eq. (39), while the last panel displays the negative 
of the source term, — (d\np/d\nT) p 5T, defined in Eq. (46). 
The response of the crust is qualitatively different in these 
cases, as we discuss below. 
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Figure 12. Same as Figure 10, but for a deep capture layer, 
Q = 95 MeV. 

For concreteness about signs, let us focus attention on 
an angular location in the star where Re{Y22{0, 4>)} is pos- 
itive, and consider a perturbation with positive 6T(r) at 
the capture layer. (Figures 10-12 are drawn to reflect this 
choice.) We imagine drawing, at this (0, <f>), a cylindrical tube 
that extends through the crust and into the fluid on both 
sides. Since we take ST to be positive at the capture layer, 
in our tube the e~ captures occur at lower-than-average Ef, 
which tends to make the capture layer relatively "underpres- 
sured" (i.e., the source term (dlnp/dT) p 8T is negative, see 
the bottom panels of Figures 10 - 12). 

In Figures. 10 - 12 the location of the capture layer, n ay , 
is identified by a sharp "kink" in £ r (for the deep capture 
layer shown in Figure 12, the kink is quite small, but its loca- 



Figure 13. Response of the crust to a smooth composition vari- 
ation, A/u e //ie = 5 X 10~ 3 , throughout the crust. Top two panels 
show vertical and horizontal Lagrangian displacements, £ r and 
£x> m meters, the two middle panels show the stresses, Ar rr and 
Ar r x, in 10 29 crgcm — 3 . The bottom panel shows the negative of 
the source term, — (dlnp/dfj, c ) p Afj, e (see Eq. [45]). 

tion is coincident with the extremum of the source term). If 
the capture layer were infinitely thin, then there would be a 
discontinuity in £ r at that location. This is because we allow 
changes of composition (i.e., relabeling of fluid elements) at 
the capture boundary, and so £r(n ay ) does not indicate the 
perturbed location of the capture layer. 

The kink in £ r is easy to understand: d^ r /dr is large 
and negative there because the crustal matter is vertically 
compressed around the capture boundary. We can estimate 
the jump in £ r at the capture layer by integrating Eq. (43a) 
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in a thin region around r lay , where the delta- function- like 
source term AS is dominant and balances dzi/d\nr. With 
this approximation, we get 



A£ r 



lay 



rh Apk B ST rflnp = A ^ Ap 
p Q d In Ef p 



(53) 



where h is the scale height at the capture layer with thresh- 
old energy Q, and Ap/p is the density jump across the cap- 
ture layer. 

Compression of the capture boundary occurs in all so- 
lutions, regardless of the capture threshold Q. Most other 
qualitative features of the solution depend on whether the 
capture boundary occurs at relatively low density, Q <C 
25 MeV, or high density, Q > 40 MeV. 

Consider first the case of a deep capture layer, shown 
in Figure 12. This is the most interesting case, since it 
generates the largest Q22- To understand the solution, we 
find it helpful to picture the displacements as built up in 
a series of "steps." (Of course, the implied time-ordering 
is fictitious; the "steps" are simply an imaginative device 
to guide intuition.) Due to the "extra" e~ captures in our 
imaginary cylindrical tube, the capture layer starts out "un- 
derpressured." Step I is vertical adjustment of the crust 
to compress the capture layer and de-compress the crust 
(d(r 2 £ r )/dr > 0) above and below it. In effect, the crust ad- 
justs vertically to "share" the "underpressure." After step I, 
the entire crust in the tube is relatively underpressured. 
Hence, in step II, the crustal matter gets "pushed in from 
the sides" at all depths. One can see this behavior in our 
solution: in Figure 12, fx is positive everywhere, so the di- 
vergence of £±(r)V a Y22 is negative (see Eq. [41]). This con- 
vergence of solid matter in from the sides causes the crust 
to be "thicker than average" in the tube; the bottom of 
the crust is pushed down and the top pushed up. Finally, 
step III: the bulging out of the crustal boundaries displaces 
the fluid in the ocean and the core, pushing it out of the 
tube. 

In summary, for deep captures, solid (elastic) matter 
is pushed into the tube from the sides, while liquid in the 
ocean and the core is pushed out. The center of mass of 
the tube also sinks somewhat. Moreover, it turns out that 
more fluid goes out of the tube than solid comes in. Both of 
these effects - the net loss of mass from the tube and the 
fact that the tube's center of mass sinks - contribute to Q22 
with a negative sign. Therefore, not only is the resulting Q22 
smaller than the fiducial estimate, Qua, but it has the op- 
posite sign! (Recall that in the picture behind the fiducial 
estimate, matter simply comes in from the sides in response 
to the "underpressure" of the capture layer. But the actual 
response of the crust is much more complicated, and numer- 
ical solutions were crucial to reforming our intuition.) 

In our deep-capture solution, £ r is of order ~ 10 m 
and fx ~ 100 m (for accretion rate of 0.5MEdd), so the 
displacements are rather large. This is clearly due to the 
smallness of the crust's shear modulus relative to pressure. 
Note also that the shear stress term 5r r ± is sizeable over 
a region roughly ~ 1 km thick. Thus, despite the fact that 
the capture boundaries are rather thin (several meters), the 
induced stresses are 'shared' by a sizeable fraction of the 
crust. This is important for the consistency of our approach. 
If we had found that all the shear stresses were concentrated 



near the capture layer, then they would certainly be large 
enough there to crack the crust, invalidating our solution. 
Even given that the stresses are shared by a large fraction 
of the crust, it is still a serious question whether the NS can 
sustain them; we address this question in detail in § 6. 

We consider next solutions sourced by a uniform 
An e /p e , since they are extremely similar to deep-capture so- 
lutions, as seen evident from comparing Figures 12 and 13. 
Again, to fix signs, take A/i c //i e to be positive inside our 
imaginary tube. The perturbation corresponds to the crust 
inside the tube being "a little underpressured everywhere." 
Compare this to the deep-capture layer case, where the crust 
starts out "very underpressured in a thin layer", but, in 
step I (described above), adjusts vertically to "share the un- 
derpressure." That is, in the deep capture layer case the 
crust adjusts vertically to resemble the uniform Ap e / fx e 
case. Of course, one can also imagine building up a smooth 
source term from a sum of delta-function-like sources. Since 
the deep capture solutions all have a similar profile, and 
since our equations are linear, the solution sourced by a 
uniform A/i c //i c must share that profile. 

We see from Figure 10 that the behavior for shallow 
capture layers is qualitatively quite different. In that case 
the crust above the capture layer sinks and compresses (the 
opposite of what we find for deep capture layers). The crust 
below the capture layer also sinks and compresses, but less 
so. Crust and fluid enter the tube above the capture layer 
and get pushed out of the tube below the capture layer. We 
find that the net effect is a positive contribution to Q22, 
though a much smaller one than the fiducial estimate. 

We see that for fixed ST, the contribution of a capture 
layer to Q22 changes sign as one goes from shallow to deep 
layers. There is an intermediate regime, Q ~ 25 — 35 MeV, 
where the displacements produce almost no contribution to 
the star's quadrupole moment. This does not mean that the 
effects from shallow capture layers in the neutron star will 
mostly cancel the effects from the deep ones; contributions 
to Q22 from deep capture layers are roughly two orders of 
magnitude larger than those from shallow layers, so the deep 
ones dominate. A complete description of how Q22 varies 
with the capture depth and other variables is presented in 
§ 5. 



5 GENERAL PROPERTIES OF 

QUADRUPOLES FROM CAPTURE 
BOUNDARIES 

In this section, we show our results for the NS quadrupole 
moment Q22, and explain how Q22 scales with physical pa- 
rameters, such as M. Our results basically validate the orig- 
inal mechanism put forward by Bildsten (1998a): small tem- 
perature perturbations in the deep crust can easily generate 
the Q22 needed for GWs to balance accretion torque. We 
defer to § 6 the question of whether the crust can actually 
sustain such large stresses without cracking or yielding. 



5.1 The Dependence of Q22 on the Neutron Star 
Parameters 

Figure 14 shows the Q22 resulting from a single capture layer 
(at e~ capture threshold energy Q). The thermal perturba- 
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Figure 14. Quadrupole values for a single capture layer, sourced 
by the temperature perturbations computed as described in § 3. 
The horizontal axis gives the threshold energy of the reaction. 
Filled squares correspond to model S (superfluid core), while open 
diamonds correspond to model N (normal core). The accretion 
rate is 0.5MEdd in all cases. The absolute value of Q22 is plotted, 
however, the sign of it is positive where the lines are solid, and 
negative where they are dotted. 



tion is due to a local nuclear heating rate that is assumed 
to vary laterally by 10% (i.e., / nuc = 0.1). Filled squares 
are for the background model with a superfluid core (gap 
energy A = 1 MeV) and open diamonds for the model with 
a normal core. The two models have the same background 
hydrostatic structure and accretion rate (M — 0.5MEdd). 

Clearly, deep capture layers generate a much larger Q22 
than shallow ones. We understand this as follows. If we ig- 
nore the vertical readjustment of the crust (as was done by 
Bildsten 1998a), then a local increase in temperature, ST, 
increases the local e~ capture rate and causes the location of 
the capture layer to shift upward by a distance Azd, given by 
Eq. (13) (see also Figure 3 and Eq. (A18) in Appendix A). In 
particular, in the outer crust, where the pressure is supplied 
by degenerate relativistic electrons, 



outer crust 



ksTST ( ST \ ( 2 

T ^^-= 30cm llO^Kj U 



,(54) 



where T ~ 10 — 20 is introduced in Eq. (12) and defined 
precisely in Appendix A. When the elastic readjustment of 
the crust is ignored, the quadrupole is just Qfld = ApAzdi? 4 
(Bildsten 1998a and see also Figure 2 in § 1.3). In the outer 
crust, where in each capture layer only 2 electrons are con- 
sumed and no neutrons are emitted, the density jump is just 
Ap/p = Ap c /p, c = 2/Z, so 



1.3 x 10 37 g cm 2 7?e 



Go 7 k) ( 



30 MeV 



Figure 15. The "sinking penalty" for the thermal perturbations 
at 0.5.Mrjdd- The absolute value is plotted, however, the sign of it 
is positive where the lines are solid, and negative where they are 
dotted. 



In the inner crust, the capture layer is rather thick, so rather 
than talk of a "jump" in density across the boundary, we 
define the fiducial quadrupole moment as 



Qfid 



<91np 



91nT 



ST i +2 
—pr dr 



. T Mc dln^o i +2 

Az d / -pr T—pr dr 

1 P dr 



(56) 



(57) 



(55) 



where F p and r^ e are defined in Appendix A; i.e., we in- 
tegrate the Lagrangian density perturbation Ap under the 
condition of Ap = and £ r = (no elastic readjustment of 
the crust, see Eq. [35]). In going from Eq. (56) to Eq. (57) 
we used Eqs. (A10), (A14), and (A18) from Appendix A. 

In the inner crust the dependence of Azd and Qfid on 
depth is more complicated than in the outer crust because of 
the variation in number of electrons captured and neutrons 
emitted in each capture layer (see Fig. 3). Also, while Q nd 
given by Eq. (56) does increase with depth for a fixed ST, 
because of the radial structure of the temperature pertur- 
bations for fixed / nU c or / C om P , Qfld must decrease near the 
very bottom of the crust. The high thermal conductivity of 
the core forces ST « at the crust-core boundary (see Fig- 
ure 8). This explains the decrease in the induced Q22 very 
close to this boundary seen in Fig. 14. 

Figure 14 shows that a single capture layer near the 
bottom of the crust generates Q22 ~ 6x 10 37 g cm 2 (//0.1), 
where / is either / nuc or / romp . There will be more than 
one capture layer near the bottom of the crust, and so the 
full Q22 for the NS is proportionately higher. The numbers 
we quote will always be for a single capture layer, so this 
multiplicative factor must be kept in mind. 
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Figure 16. Quadrupolc values for a single capture layer, with 
Aza = 100 cm plotted as function of the threshold energy of 
the reaction. The lines marked with filled squares correspond to 
the full crust, while the open diamonds represent the case where 
the crust is melted around neutron drip. The line marked with 
upward-pointing triangles corresponds to the crust with the shear 
modulus fi artificially decreased by a factor of 2 from the value 
given by Eq. (Al), while the line marked with downward-pointing 
triangles corresponds to the case where fi is 2 times larger than 
the value from Eq. (Al). Finally, the line marked with open circles 
corresponds to the case where there is a 5% density jump at the 
crust-core boundary (i.e. pi/p s = 1-05). The absolute value of 
Q22 is plotted, however, the sign of it is positive where the lines 
are solid, and negative where they are dotted. 

Though the scaling of the fiducial estimate, Qfld, is help- 
ful for our understanding, it neglects an essential piece of 
physics. As described in § 4.3, the crust prefers to sink in re- 
sponse to the shift in capture layers, and this reduces the ac- 
tual quadrupole moment significantly below the fiducial esti- 
mate. In Figure 15 we plot this "sinking penalty", (522/Qfld, 
for the same model as in Figure 14. As one can see, the 
penalty is quite large (~ 20 — 50) for the shallow capture 
layers (Q < 40 MeV), but is only ~ 5 — 10 for the capture 
layers in the deep crust. The fact that the sinking penalty de- 
creases with greater depth, while Qfld increases, means that 
deep capture layers are the dominant contributors to Q22- 
We find that transverse temperature contrasts of 10 6 — 10 7 K 
are sufficient to generate Q22 ~ 10 37 — 10 38 g cm 2 . This is 
much smaller than the ~ 10 8 K contrast originally required 
by Bildsten (1998a) for the shallow capture layers. 

In Figure 16 we explore the effects of the magnitude 
of the shear modulus fj,, the physical extent of the crust, 
and the possible density discontinuity at the crust-core in- 
terface. In order to simplify the discussion, we computed 
the quadrupole moments resulting from a fixed shift Az d = 



100 cm of the capture layers, rather than from a temperature 
perturbation ST that has a non-trivial radial dependence 
(see Appendix A for a detailed discussion). Filled squares 
show Q22 for the standard model used throughout this sec- 
tion. The Q22 values are different from those in Figure 14 
since in Fig. 16 the perturbations are normalized to yield 
the same vertical shift Azd of the capture layers regardless 
of its position in the star. 

As shown by Brown (2000), at high M the energy in- 
put from crustal nuclear reactions can melt the crust near 
neutron drip. In this case there will be a liquid layer in the 
middle of the crust, and our outer elastic boundary condi- 
tion needs to be applied there. The line marked with open 
diamonds in Fig. 16 shows Q22 as function of capture layer 
depth for this case. The quadrupole moment due to deep 
capture layers is virtually unaffected (differs by ~ 10%) by 
the presence of a melted layer near neutron drip. 

The lines marked with triangles survey the dependence 
of Q22 on the crust shear modulus p,. Doubling fi com- 
pared to the canonical value roughly doubles Q22 (the line 
marked with downward-pointing triangles), and halving p 
(line marked with upward-pointing triangles) roughly halves 
£?22. This is consistent with the general dependence of the 
quadrupole moment on the shear modulus, as derived in § 6. 

Throughout this paper we have assumed that the den- 
sity is continuous across the crust-core boundary. However, 
this is not certain (see Pethick & Ravenhall 1995 for a re- 
view), as many models for the transition contain density 
jumps. Though gravitationally stable, this density discon- 
tinuity adds to the restoring force when the crust pushes 
down into the core. The line marked with open circles shows 
Q22 for a NS model where the liquid just inside the core 
is 5% denser than the crust immediately above it (i.e., 
pi/ ps = 1.05). In this case, the rippled surface of the crust- 
core boundary provides some extra restoring force that re- 
duces Q22- 

5.2 Dependence of Q22 on Accretion Rate, and 
the Overall Picture 

We now turn our attention to the scaling of Q22 with M. 
This behavior is important for comparing our results to the 
observed distribution of NS spin frequencies in LMXB's, 
since in our picture the close similarity in spin frequencies 
(all ~ 300 Hz) over a large range in M implies that Q22 
scales (roughly) like M 1 ^ 2 . In Figure 17, we show the Q22 
induced by 10% perturbation in nuclear heating (/ nu c = 0.1, 
lines marked with filled triangles) and by 10% perturbations 
in Z 2 1 A (/comp = 0.1, lines marked with open circles), as 
a function of M. Solid lines correspond to models with a 
superfluid core, while dashed lines indicate models with a 
normal core. The dotted line displays, as a function of M, 
the quadrupole required for spin equilibrium at v a = 300 Hz 
(see Eq. [4]). 

It should be clear from this figure that the Q22 gener- 
ated by a single capture layer can easily account for the spins 
of the accreting neutron stars. Given that there are several 
capture layers in the deep crust, the required asymmetry 
in either the heat sources or compositions is thus propor- 
tionately lower than 10%. Hence, despite the reduction in 
Q22 (below Qfld) due to the vertical readjustment of the 
crust, the basic mechanism of electron capture induced den- 
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Figure 17. The quadrupole moment Q22 due to a single cap- 
ture layer with threshold energy Q = 95 MeV. Solid lines de- 
note the results for the model with a supcrfluid core (gap energy 
A = 1 MeV), while dashed lines denote the values for the model 
with a normal core (A = 0). Lines marked with open circles are 
for the quadrupole moment sourced by the composition perturba- 
tions, / comp = 0.1, while lines marked with filled triangles denote 
the quadrupoles sourced by the variation in the local heating of 
the crust, / nuc = 0.1. Finally, the dotted line is the relation given 
by Eq. (4), i.e., the quadrupole moment necessary for spin equi- 
librium at u s = 300 Hz as a function of M. 



sity jumps can indeed account for the spin rates. It is also 
evident from Figure 17 that the case of superfluid core with 
ST sourced by lateral composition variations / C om P (solid line 
marked with circles) can be eliminated as a viable source of 
the quadrupole moment, especially at low accretion rates, 
as /comp would need to be close to unity. 

Remarkably, for accretion rates M < 0.2M E dd, Q22 for 
a constant / nuc or f CO mp scales as M 1 ^ 2 , i.e., in the same way 
as the required quadrupole moment, Eq. (4). At high accre- 
tion rates, the temperature profile in the inner crust is set 
by the balance of the heat input due to compression-induced 
nuclear reactions near neutron drip and neutrino emission 
from the crust and the core (Brown 2000). On the other 
hand, at low accretion rates, the effect of heating near neu- 
tron drip and crustal bremsstrahlung is diminished, which 
leads to a more direct dependence of the temperature profile 
on the accretion rate and the outer boundary temperature. 

As we said before, the NS's total Q22 depends on the 
magnitude of the asymmetry (i.e., / nuc or / C om P ) and the 
number of capture layers, i.e., the curves showing Q22 for 
different models in Figure 17 can be shifted up or down by 
overall factors. Despite this uncertainty in the prefactor of 
<?22(M), our calculation has no uncertain parameters that 



Figure 18. The amplitude of flux perturbations SF r at the top 
of the crust that corresponds to a Q22 due to a single capture 
layer given by equation [4], i.e. a quadrupole moment that ensures 
equilibrium v s = 300 Hz at a given accretion rate. The legend is 
same as in Figure 17. &F r is proportional to the number of capture 
layers contributing to Q22, so with many capture layers in the 
deep crust the required 8F r is reduced accordingly. According to 
Eq. (4), SF r oc v s , i.e., for equilibrium spin of 600 Hz the 
estimated 5F r is reduced by a factor of 5. 

would alter the calculated scaling of Q22 with M for an M- 
independent / nuc or / comp . 

As discussed in § 3.3, the thermal perturbation calcula- 
tions also predict a modulation in the X-ray flux exiting the 
NS crust. The temperature perturbations in the deep crust 
that displace the capture layers also give rise to a lateral 
variation of the flux at the neutron star surface. These hot 
and cold spots, when moving in and out of the view of the 
observer due to the rotation of the neutron star, generate a 
modulation in the persistent emission. 

The perturbed flux 5F r arising from either / nuc or / CO mp 
equal to 10% was shown in Figure 9 of § 3.3. However, as 
we see in Figure 17, a 10% / nuc or / comp perturbation in a 
star with a normal core produces Q22 from a single capture 
layer that is too large at low M and too small at high M. If 
the equilibrium spin frequency is exactly independent of M, 
then for our deformed capture layer model to accord with 
these observations, / nuc or / comp must vary with M. 

In Figure 18 we show 5F r , normalized by the accre- 
tion flux Fa ~ 200 Me V/4itR 2 /accreted baryon, under the 
assumption that the Q22 due to a single capture layer at 
Q = 95 MeV is exactly that needed to set v s = 300 Hz 
(Eq. [4]). In other words, we adjust / nuc and / CO mp as a func- 
tion of M to get the required Q22- The actual flux asymme- 
try in the NS is of course lower than that shown in Figure 18 
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by the factor of (roughly) the number of deep capture layers, 
as additional capture layers in the deep crust will increase 
the quadrupole moment for the same 5T/T. 

The ratio SF t /Fa shown in Fig. 18 presumes that accre- 
tion is steady. However, there are systems, such as neutron 
star transients, that accrete only episodically, with time- 
averaged accretion rates (M) < 10 _11 M Q yr _1 . When ac- 
cretion halts in these systems, the thermal emission from 
the surface should be visible directly (Brown et al. 1998). 
Hence, in quiescence, the modulation is 5F r / F(r top ) ~ 
1(T 2 (//10%), rather than SF r /F A . 

Therefore, our thermal calculations lead us to predict 
a certain level of modulation in the X-ray luminosity of 
LMXB's. Can this modulation be detected? Possibly. We 
have not considered how the ocean would respond to such 
lateral flux asymmetry emerging from the crust. One could 
imagine transverse flows being generated from the resulting 
transverse temperature gradients. Presuming that the flux 
at the NS photosphere maintains the asymmetry of the flux 
at the top of the crust, it can provide a critical window onto 
processes in the deep crust. The bright source Sco X-l would 
be a good place to start, as it has the highest count rate. 
The best current limit on a coherent pulse from Sco X-l is 
< 1% (Vaughan et al. 1994). Our capture-layer mechanism 
for quadrupole moment generation might then be tested by 
observations that are accessible today, without having to 
wait 5 — 10 years for gravitational wave detectors to reach 
the requisite sensitivity. In addition to the excitement of 
finding a coherent pulse in the X-rays, such discovery would 
substantially reduce the parameter space that would need to 
be searched for the gravitational wave detection. 



6 SHEAR STRENGTH OF THE CRUST AND 
THE MAXIMUM POSSIBLE CRUSTAL 
QUADRUPOLE 

In the previous sections we described how lateral pressure 
gradients due to composition asymmetries can deform a neu- 
tron star crust and computed the resulting quadrupole mo- 
ments. However, there is a maximum degree of deformation 
that a crust can sustain, which is set by its yield or break- 
ing strain <7 max . The yield strain places a fundamental limit 
on Q22 of a neutron star. Unfortunately, yield strains of 
even terrestrial materials, let alone neutron star matter, are 
poorly understood. In § 6.1 we summarize the current folk- 
lore regarding o- max . 

Wavy electron capture layers are a very special way of 
straining the crust. However, absent strong magnetic fields, 
a (static) neutron star's quadrupole moment must be gen- 
erated by shear stresses in the crust, no matter how the 
stresses arise. Imagine somehow turning off the crust's shear 
modulus. Then the matter equations are simply the Euler 
equations for a perfect fluid, and then Q22 must be zero in 
static equilibrium. These considerations suggest that there 
should be an expression for Q22 that involves the shear-stress 
forces only. In § 6.2 we derive such an expression, Eq. (64). 

Eq. (64) allows us to estimate the quadrupole moment 
that results from a "typical" strain amplitude. It also gives 
us an upper limit on Q22, set by the yield stress (derived 
in § 6.3). Given the tremendous uncertainty in a ma x, this 
formula allows us to lump all our ignorance into a single 



parameter. Note that our expression (64) for Q22, and the 
resulting upper limit (69), do not depend on several other 
assumptions made in the rest of this paper: that the "non- 
stressed" state of the crust is spherical, that the stress-strain 
relation is in the linear regime, or that the perturbation 
is sourced by temperature or composition gradients. It ap- 
plies regardless of the crust's detailed evolution, which could 
have undergone an arbitrary amount of creep or cracking. Of 
course, our formula does not apply when Q22 is generated 
by some force other than shear stresses, such as magnetic 
fields. 

In § 6.4, we use the formalism developed earlier to eval- 
uate the strain induced in the crust by the wavy e _ capture 
layers or uniform composition gradients. We compare the 
strains to the crust's maximum strain level and discuss the 
implications. 

Finally, since the analysis in this paper ignores the self- 
gravity of the perturbations, in § 6.5 we estimate the size of 
the correction when self-gravity is included. 

6.1 Maximum Strain Level for a Neutron Star 
Crust 

Throughout this paper we have made a simplifying assump- 
tion that the response of the crust to the density and pres- 
sure perturbations is purely elastic. However, the response 
of solids to applied stresses is more complicated than that. 
There are two related issues here. First, solid materials be- 
have elastically only up to some maximum strain a- max , be- 
yond which they usually either crack or deform plastically. 
Upon relieving the stress, the solid does not return to its 
initial shape. Second, even at strains well below the yield 
strain, a solid that has been strained for a very long time 
tends to "forget" its former equilibrium shape; i.e., the equi- 
librium shape undergoes irreversible relaxation or "creep." 
This behavior is called viscoelastic: solids respond elastically 
on short timescales, while on very long time scales they be- 
have more like very viscous fluids. As argued in § 4, we ex- 
pect that viscoelastic relaxation, while perhaps important 
to the detailed picture of "NS mountain-building," will not 
drastically change our conclusions about the likely magni- 
tudes of NS crustal deformations. However, the question of 
the maximum strain that a neutron star crust can sustain is 
a crucial one for our work, as is the question of what happens 
when this strain is exceeded. 

Yield strains a" max of even ordinary materials, let alone 
neutron star crusts, are very hard to predict theoretically. 
For perfect one-component crystals, simple theoretical con- 
siderations lead to <7max ~ 1CP 2 — 1CF 1 (Kittel 1956). How- 
ever the maximum strain of most solids is determined by the 
motion of dislocations and other defects. Early discussions 
by Smoluchowski (1970) placed <7 max of neutron star crusts 
in the range of 10~ 2 to 10~ 5 by analogy with a variety of 
terrestrial materials with chemical and lattice imperfections. 
Ruderman (1991) puts the maximum strain in the range of 
10~ 4 — 5 x 10~ 3 , which implies a large number of cracking 
events in the lifetime of a spinning down young pulsar. 

Recent breakthroughs in the understanding of the soft- 
gamma ray repeaters as magnetars may shed some light on 
our problem as well. Thompson & Duncan (1995) argued 
that the energy release that powers these events is mag- 
netic in nature, but that the events are triggered by "crack- 
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ing" of the crust, which accumulates strain as the magnetic 
field evolves. The maximum energy of such events implies 
(in a model-dependent way) that the maximum strain is 
in the 1CP 3 range. Phenomenological support for the idea 
that crust cracking is the origin of the bursts comes from 
the rather striking similarity of the power-law distribution 
of burst energy with that of terrestrial earthquakes (Cheng 
et al. 1996). 

However, cracking is just one way in which the crust can 
relieve the applied stress. Smoluchowski & Welch (1970) ar- 
gued that hot pulsar crusts might undergo large amounts of 
plastic or ductile deformation and not crack nearly as of- 
ten, an effect more recently argued by Link et al. (1998) as 
arising from high pressures (namely, p ^> ft). Indeed, terres- 
trial materials undergo a brittle to ductile transition when 
pressures become much greater than their shear modulus 
(Turcotte & Schubert 1982). Thus, there is a large range of 
options, none of which we are confident to exclude. 



6.2 Integral Expression for Q22 in Terms of Shear 
Stresses 



We write the stress-energy tensor of the solid as 

Tab — —PQab + tab , 



(58) 



where t a t is the (trace-free) shear stress tensor of the crust. 
Consider the "deformed" star as a spherical star plus a small 
perturbation, and consider t a b to be a first-order quantity, 
so that 

$T ab = -8pg ab + t ab . (59) 

Equilibrium between gravity and hydro-elastic forces implies 

V a 5T ab = 8pg{r)f b , (60) 

where we have neglected the influence of the perturbation on 
the star's gravitational field (the Cowling approximation). 
We expand the perturbation in spherical harmonics: 

tab = t rr (r)Yi m (f a fb - ie a6 ) (61) 

+ t r ±(r)f a b + tA(r)(A ab + ^Yi m e a b), 

where e a b, fab and A a b were defined in Eq. (40). The above 
expansion is automatically trace-free. We have left out 
terms in t ab proportional to {f a e bcd + f b e acd ) r d V c Yi m or 
(V Yi m e&cd + VbYi m e acd ) r d V c Yi m because, having oppo- 
site parity, they decouple from the other shear stress terms 
and cannot generate a quadrupole moment. 
Projecting Eq. (60) along f , we obtain: 



g(r) 



dSp dt rr 3 _ f3 
dr dr r r 



(62) 



We now replace the Sp term on the right side of Eq. (62) 
in favor of shear stress terms by projecting Eq. (60) along 

2): 



V Vim, which yields (specializing to I 
Sp(r) 



1 

~2 U 



3 



1 r dt r ±_ 
3 tA + P^ 



(63) 



Using Q22 = J 8pr 4 dr, and integrating by parts (using the 
fact that the shear stress vanishes above and below the crust 



(see § 4.2), we obtain: 

Q22 = - U) Ur + i (6 - U) t A (64) 

dr, 



8 



317+ \ll 2 



3 r ~dr~ 



where U is defined in Eq. (44). This expression gives the 
quadrupole moment of the crust, so long as it is in hydro- 
elastic balance. Since this expression involves only shear 
stresses, it also makes clear that (static) perfect fluid stars 
cannot have a quadrupole moment. 



6.3 Maximum Q22 Set by Crustal Yield Strain 

We now reinterpret Eq. (64) in terms of strains, rather than 
stresses, by defining a ab = t a b/fi, where ft is the shear mod- 
ulus. We define a by a 2 — ^a a bcr ab , and we assume that the 
crust will yield when a > a" max . This criterion for yielding is 
called the von Mises criterion (Turcotte & Schubert 1982). 
Different empirical criteria are sometimes adopted (such as 
the Tresca criterion, which depends on the difference be- 
tween the maximum and minimum eigenvalues of o a b), but 
for our purpose they would all give similar answers, and the 
von Mises criterion is easier to use. 

Note that U = 4ivr 3 p/Mr < 1 in the crust, so all the 
coefficients of the stress components in Eq. (64) are posi- 
tive. Therefore, we can compute the upper bound on Q22 by 
making the stress components t rr , t r ±, and tA as large as 
possible, subject to a a bO~ ab < 2a^ nax . Of course, the crustal 
yield strain could vary with density, but given that its value 
is so uncertain, we shall simply take it to be some constant, 
characteristic of the entire crust. 

Using our expansion (61) of t a b , we write^ 

(Jab<J ah = ^a 2 r [Re(Yi m )] 2 + a 2 ± [Re(f a b)] 2 (65) 
+ cr 2 K [Re(k ab + ^Yi m e ab )] 2 . 

For I — m = 2, the following identity holds among our basis 
tensors: 

3 r „ „, . 3 r „ ,. „ 9 . 9 r „ , 1„ N1 o 15 



'-[Re{Y lm )] 2 + -[Re{fab)] 2 + l\Re{Kb+\Y lm eab)] 2 



32tt 



(66) 



Therefore a will attain its maximum value, a max , at every 
point in the crust — i.e., the crust will be everywhere stressed 
to the maximum — if we set 



/32tt\ 1 /2_ 
<j rr — \ 15 / ^ max ' 


(67a) 


/167ry/ 2 

C"r± = 1 — g— 1 Cmax, 


(67b) 


/m-KX 1 ' 2 

C"A = 1 — — J Crnax- 


(67c) 



With this substitution, all the terms in the integrand in 
Eq. (64) share a common prefactor, fir /g. Moreover, since 
U <C 1, the integral in Eq. (64) can be expressed as Q22 = 
l^max-/", where the / is the integral 



M Recall that we defined all our variables to be real. 
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-I 



pr 3 dr 



(68) 



which depends strongly on M, R, and the location of 
the crust-core boundary, while 7 is a numerical prefactor 
that depends very weakly on those parameters. / is ap- 
proximately given by / ~ (2/7) (p/p) (pb/g)R 3 AR, where 
Pb is the pressure at the crust-core interface, AR is the 
crust thickness, and (p,/p) is a (suitably) weighted av- 
erage of the shear modulus in the crust. At densities 
well above neutron drip, the shear modulus is p/p ~ 
10~ 3 (Z/20) 2 (88/,4) 4/3 p~ 4 1/3 , assuming the neutron fraction 
X n = 0.8. The pressure is approximately ideal degenerate 
neutron pressure, p n oc p 5 ^ 3 , and AR is roughly 5/2 times 
scale height at pb, i.e., AR ~ (5/2)pb/ pbg- Hence, we ex- 
pect the integral I to scale as / oc p 2 R 7 M~ 2 . By fitting to 
our detailed numerical calculations, we get the more precise 
scaling I = 3.61 x 10 37 g cm 2 p\f 
terms, we find 



' M^l' 2 R$ ' Aa . Collecting 



Q max = 1.2 x 10 38 g cm 2 (|p 

/ 1-X n \ 4 / 3 /x n \ 6 / 3 
X I 0.2 J \0.8J 



MJ : 2 V20J I 



88 \ 

a) 



Pb 



2.1 x 10 14 g cm- 



4/3 



(69) 



where the values of Z, A, and X n are understood to be 
averaged over the crust, with a weighting prefactor pr 3 /g, 
which is heavily biased to select the values near the bottom 
of the crust. In Eq. (69) the prefactor and the scalings for pb, 
M, and R are from numerical calculations, while the terms 
in Z, A, and X n simply come from our analytical scaling 
formula for the shear modulus. 

Eq. (69) is very powerful. It shows that, so long 
as only elastic forces are important, the maximum 
quadrupole moment attainable for a NS crust is Q22 ~ 
10 38 g cm 2 (<T max /10 -2 ), no matter how the strains arise. The 
exact upper limit does depend on the NS radius, however. 
Note that for a given NS EOS and crust composition, the 
upper bound on Q22 is smaller for heavier NS's (which have 
smaller radii) . We also emphasize that even to approach this 
upper limit, almost all the strain must be in the Y22 spher- 
ical harmonic; strain in other harmonics pushes the crust 
closer to the yield point without contributing to Q22- 

Also note that, besides providing an upper limit, our 
formula (64) provides an estimate of the Q22 that can 
result for a given level of strain in the crust: Q22 ~ 
10' g cm ((0-22} /10" 2 ), where {a 22 ) is some (appropri- 
ately) weighted average of the "22-piece" of the crustal 
strain. 

Note that our formula, suitably interpreted, also holds if 
the crust's stress-strain relation is in the non-linear regime. 
In that case, we just define an "effective" strain u^f by 
a lb = tab/p, where p is the shear modulus valid in the 
linear regime. Then Eq. (69) continues to hold if we replace 
by (J max, the maximum value of the effective strain. 



6.4 Crustal Strain Due to Electron Captures 

What strains are induced by wavy e~ capture layers or a 
uniform composition gradient? Expanding the shear tensor 
as functions of r times (same-parity) basis tensors built from 
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Figure 19. Angular maximum values of the shear strain com- 
ponents max(9 0-){iT a i,} as a function of the position in the crust 
due to the motion of a deep capture layer (Q = 95 McV) in re- 
sponse to an / nuc = 0.1 thermal perturbation. The background 
model has a normal core, the accretion rate is 0.5 Mgcld! an d the 
quadrupole moment is Q22 = 8.8 X 10 37 g cm 2 . Solid line de- 
notes (15/327r) 1//2 oy r (r), dashed line indicates (b/SnY^Or^r), 
and dash-dotted line shows (5/487r) 1 / 2 o'A( r )) where a a b(r) are as 
defined in Eq. (70). 



Yi m , as in Eq. (38), and then expressing the radial functions 
in terms of our variables zt, we have 



a rr (r) 
a r± (r) 



4 fdjr <ir\ | 2/Sei 

3 V dr r J 3 r 



4 dzi 2 2 

7T "Ti 1" Z3 

3 a m r 3 



d 5x 

r ^ h 

or r r 

2 ^L = 2f 3 2 z 3 . 
r 



P 2 z 4 



(70a) 

(70b) 
(70c) 



The rr-component of strain at any point in the crust is 
then a rr (r,d,4>) = a rr (r)Re{Yi m (6, 0)}, and similarly for 
the other components. In particular, for I = m = 2, the 
maximum value of the angular factor is (15/327r) 1/ ^ 2 = 0.39 
for a rr , (5/8tt) 1/2 = 0.45 for a r± , and (5/48tt) 1/2 = 0.18 
for cta- In other words, max{aj L [Re(A ab + (l/2)Yi m e ab )] 2 } = 
5/48-7T, and similarly for a rr and a r ±. 

Figures 19 and 20 show the components of shear strain 
induced by the wavy capture layers and smooth composition 
gradients, respectively. In Figure 19 the temperature pertur- 
bation comes from a thermal model with M = 0.5AfEdd, and 
/nuc = 0.1; the resulting temperature perturbation at the 
Q = 95 MeV capture layer is 8T = 9 x 10 6 K (8T/T m 2%, 
see Figure 8) , and the deformation of this layer leads to 
Q22 = 8.8 x 10 37 g cm 2 . The magnitude of the composition 
gradient in Figure 20 is Ap e /p c = 5 x 10~ 3 , which results 
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Figure 20. Same as Figure 19, but for displacements sourced by 
a smooth composition gradient with A/^ c /Vc- = 5 X 10~ 3 . The 
quadrupole moment in this case is Q22 = 1.2 X 10 38 g cm 2 . 



in Q22 = 1.2 x 10 38 g cm 2 . Note that for these fiducial per- 
turbations, the resulting Q22 values are basically what is 
needed for gravitational waves to balance accretion torque 
in LMXB's, with M « MEdd- Of course both the strains and 
the quadrupole moment scale linearly with / nuc or A/j, e //j, e . 

From these figures it is evident that a rr is much big- 
ger than a r ± and a a- Our mechanism is thus not optimally 
efficient at producing Q22, since the crust ends up being 
much more strained at the equator than at the poles. (I.e., 
a solution that is near the yield strain at the equator is still 
relatively unstrained near the poles, and so the polar regions 
could be "doing more" to hold up the quadrupole.) In fact, if 
a rr is the dominant stress component, then Eq. (69) becomes 
Q22 = 3 x 10 37 g cm 2 ((<r rr ) /1CF 2 ). The strain oy r near the 
bottom of the crust in these two models is ~ 4 — 5 x 1CF 2 , 
which shows that this estimate accords nicely with our nu- 
merical solutions. 

The second thing to note about Figures 19 and 20 is 
that a rr > 1CF 1 near the top of the crust. Thus our lin- 
ear elastic model is inconsistent there, in the sense that our 
solutions must certainly exceed the yield strain. To us, this 
suggests that our assumed level of composition inhomogene- 
ity would drive a continual plastic deformation of the crust 
there. Although we are therefore using the "wrong physics" 
to describe the top of the crust (i.e., we should be using 
equations that somehow incorporate plastic deformation), 
we believe this will not lead to significant errors in our main 
results-the values of Q22- This is because, for the solutions 
shown in Figures 19 and 20, the top of the crust contributes 
only a small fraction to the total Q22- For example, if we 
repeat these calculations, but with top of the crust set at 
around neutron drip (r = 9.2 x 10 5 cm), we find that Q22 dif- 



fers by ~ 10% from the value obtained with a full crust (see 
Figure 16). As long as most of the crustal mass is strained 
below the yield level, our results for Q22 are reasonably ac- 
curate. 



The fact that a smooth composition gradient induces a 
Q22 = 1.2 x 10 38 g cm 2 ((A// e /Me)/0.5%) has interesting as- 
trophysical implications. As is evident from Figure 20, even 
such small composition gradients induce sizeable strains in 
the crust. Moreover, unlike the quadrupole moments pro- 
duced by deformed capture layers, which would be wiped 
out in the time it takes to accrete the mass in a capture layer 
(Bildsten 1998a), if the thermal gradient is turned off, the 
uniform composition gradient can only be eliminated by re- 
placing the entire crust; i.e., accreting ~ O.O5M0 of material. 
Hence, even when accretion ceases, the crust of an LMXB 
neutron star is likely to have a remnant quadrupole moment. 
In transiently accreting systems, this remnant quadrupole 
moment may set an upper limit on the spin frequency. 

In general, it is clear that at high accretion rates (;> 

0.5MEdd), the quadrupole moment needed to balance ac- 
cretion requires the crustal strain to be ;> 10~ 2 . This is 
probably higher than the yield strain, so if such equilib- 
rium prevails, it seems the entire crust must be in a state of 
continual plastic flow. Assuming that accretion continually 
deforms the entire crust by the above mechanism, then the 
stresses are likely to stay near the yield value. This can pro- 
vide a natural explanation for the similarity of spin frequen- 
cies in near-Eddington accretion rate systems. We have not 
attempted to model the resulting plastic flow, but that may 
be worth pursuing. We caution, however, that this picture 
of inhomogeneity-driven plastic flow is based solely on the 
folklore regarding yield strains of materials, as no definitive 
calculations of a max exist, and estimates are typically based 
upon extrapolating experimental results for ordinary terres- 
trial materials by > 10 orders of magnitude. Our approxi- 
mations have allowed us to quantify quite clearly how large 
the crustal stresses must be for gravitational wave emission 
to be appreciable in accreting neutron stars. How the crust 
responds to such high stresses is a problem for future re- 
search. 



Since the quadrupole moment scales linearly with the 
crustal strain, in lower accretion rate systems (< 0.5M E dd), 
the required strains are correspondingly lower (~ 10~ 3 — 
10~ 2 ). Moreover, we showed in § 5 (see Figure 17), that 
for a fixed composition asymmetry (/ nuc or / C om P ) Q22 has 
exactly the right scaling with M to balance the accretion 
torque by mass quadrupole gravitational radiation at a fixed 
spin frequency, independent of the accretion rate. On the 
other hand, for high accretion rate systems (> 0.5MEdd), 
in order to explain the spin frequency clustering at exactly 
300 Hz, our mechanism requires that the asymmetry in the 
crust (/nuc or / C om P ) correlate with M in a well-defined way 
(see Figure 17). Alternatively, if / CO mp and / nuc are the same 
as in lower accretion rate systems, then we would expect 
brighter LMXB's to have higher spin frequencies. Given the 
uncertainty in the spin frequency measurements for these 
sources, such a possibility cannot be ruled out at present. 
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6.5 Correction for Self-Gravity of the 
Perturbations 

In deriving the maximum quadrupole formula, Eq. (64) 
and (69), as well as in the rest of this paper, we have ne- 
glected the changes in the gravitational force due to the per- 
turbation itself — i.e., we have neglected the deformation's 
self-gravity (the Cowling approximation) . Very roughly, one 
would expect these self-gravity corrections to increase Q22 
by a fractional amount proportional to /crust //ns ~ 0.1. 
However, including the gravitational potential perturbation, 
5$, allows non-zero Sp in the interior of the stellar core, 
something that was forbidden within the Cowling approx- 
imation. For perturbations with low I, the self-gravity ef- 
fect turns out to be much larger than O (/crust //ns)- We 
now derive a formula for Q22 similar to Eq. (64), but in- 
cluding the self-gravity of the deformation. However, while 
Eq. (64) was exact (within the Cowling approximation) , our 
improved version will include the effects of self-gravity only 
in an approximate way. 

Including Newtonian self-gravity, our expressions (58), 
(59), and (61) for T ab , Sr a b, and t ab remain valid, but now 
the equilibrium 8r ab satisfies 



V a 5T ab = 5pg(r)f b + pV b 5$. 
Projecting Eq. (71) along f b , we obtain: 



dp = 7 7 



dSp ^ dtrr 

dr dr 



3 

+ -t r 

r 



— tr± 

r 



(71) 



£^i. (72) 
g dr 



Replacing Sp in (71) by projecting Eq. (71) along V 6 Yi m , 
and integrating by parts to eliminate radial derivatives of 
the shear terms, we obtain 



Q22 = rhs of (64) + / — ^f-S^dr 
g dr 



(73) 



where the "extra" term ' in Eq. (73) is an integral over 
the entire star, not just the crust. Deformations in the 
crust change the potential throughout the star, so now Sp is 
nonzero everywhere. 

Equation (73) is exactly true, but not very useful with- 
out knowing <$$(r). We introduce an approximation that 
allows us to obtain a closed- form expression for Q22- We 
expect that for large quadrupole moments, most of the den- 
sity perturbation lies near the bottom of the crust. As an 
approximation then, we use the <$$(r) appropriate to a thin 
deformed shell at radius rbot, the location of the crust-core 
boundary: 



**(r) = -^Q 2 2 
5 



2 /riot r < r bot 
r > rbot- 



1/r 3 



(74) 



Plugging Eq. (74) into (73) we obtain 



Q2 



where 



(4 - (j) trr + \(6-U)tA 



3 3 dr 



dr, (75) 



T = 



47T 

"~5~ 



dp 



m(r) dr 



dr + 



dp 



m(r) dr 



dr 



.(76) 



Note that T is manifestly positive (since dp/dr < 0), so the 
factor (1 — !F)~ coming from self- gravity always leads to 
an enhancement in Q22 over the value given by our formula 
(64). Typical values of T are 0.2 — 0.5, depending on the ex- 
act core model. Hence, we expect that including self-gravity 
in the full calculation will enhance the resulting Q22 by a 
factor of 25% - 200%. 



7 CONCLUSIONS 

We have investigated whether the accretion-driven spin- 
up of neutron stars in low-mass X-ray binaries (LMXBs) 
can be halted by gravitational wave emission due to 
mass quadrupole moments generated in their crusts. The 
quadrupole moment needed to reach this spin equilibrium is 
Q22 ~ 10 37 - 10 3S g cm 2 (see Eq. [4]) for the relevant accre- 
tion rates in LMXBs, 10~ 1() -2 x 10~ 8 M Q yr -1 . How to form 
and sustain a quadrupole this large is the main problem we 
addressed. We have undertaken a series of calculations that 
substantially extend the original idea of Bildsten (1998a) 
that electron capture reactions can deform the crust by large 
amounts. The major results of our work are as follows: 

(i) By self-consistently solving the elastic equilibrium 
equations (§ 4 and 5) we have found that the predomi- 
nant response of the crust to a lateral density perturbation 
is to sink, rather than move sideways. For this reason, the 
quadrupole moments due to temperature-sensitive e~ cap- 
tures in the outer crust (i.e., the case considered by Bildsten 
1998a) are actually much too small to buffer the accretion 
torque. However, a single e~ capture layer in the deep inner 
crust can easily generate an adequate mass quadrupole. This 
requires lateral temperature variations in the deep crust of 
order <^ 5%, and the realistic case of multiple electron cap- 
ture layers requires a proportionately smaller temperature 
variation. Because of the much larger mass involved in gen- 
crating Q22 in the inner crust, the temperature contrasts 
required are only ~ 10 6 — 10 7 K, rather than ~ 10 8 K orig- 
inally envisioned by Bildsten (1998a). Alternatively, a 0.5% 
lateral variation in the charge-to-mass ratio can generate a 
Q22 sufficient to balance the accretion torque even in the 
absence of a temperature gradient. 

(ii) Our thermal perturbation calculations show that the 
temperature variations required to induce a quadrupole 
moment this large can easily be maintained if the 
comprcssionally-induced nuclear reactions around neutron 
drip inject heat with about ~ 10% lateral variations. Lat- 
eral variations of the same magnitude in the average Z 2 1 A of 
nuclei in the crust can also maintain a similar temperature 
asymmetry. This is despite the strong thermal contact with 
the isothermal core of the neutron star (see § 3). However, 
if accretion halts or slows considerably, then these temper- 
ature variations will be wiped out in a thermal time, i.e., a 
few years at the crust-core interface (Eq. [31]). In this case 
the e~ capture boundary deformations will be smoothed 
out in the time it takes to accrete the mass in a capture 
layer, ~ 2.5 x 10 6 yr (M/10~ 9 M yr 1 ) near the crust-core 
boundary. 

(iii) While it is not possible to estimate the size of the 
compositional or nuclear heating asymmetries a priori, the 
temperature variations in the deep crust lead to lateral vari- 
ations in the persistent thermal flux emerging from the neu- 
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tron star. Since the neutron star is spinning, a certain Q22 
implies an amplitude of the modulation of the persistent X- 
ray flux (see Figures 9 and 18). Though a small effect, these 
periodic variations can be searched for observationally. De- 
tection of such modulation would help tremendously in the 
search for gravitational wave emission. 

(iv) If the size of the heating or composition asymmetry 
is a constant fixed fraction, then we showed that for M < 

0.2MEdd the scaling of Q22 with M is just that needed for 
all of these low accretion rate NS's to have the same spin 
frequency (see Figure 17). 

(v) Quadrupoles this large require a strain in the crust of 
order 1CP 3 — 1CP 2 (depending on the accretion rate, with 
strains exceeding 1CP 2 for M ;> 0.5MEdd, see discussion in 
§ 6.4), regardless of the detailed mechanism for generating 
the strain. 

(vi) We have derived a general relation between the max- 
imum quadrupole moment Qmax that a crust can support 
via elastic deformation and its breaking strain (see § 6). 
Our relation (Eqs. [64] and [69]) is more complete than pre- 
vious work and is widely applicable. In addition to determin- 
ing Qmax, this relation allows one to robustly estimate the 
quadrupole moment for a given level of strain, even in the 
nonlinear regime, regardless of the amount of plastic flow or 
relaxation that the crust has undergone. 

Our work has thus clarified many of the important 
outstanding questions for the hypothesis that gravitational 
wave emission due to a crustal quadrupole moment can 
buffer the accretion torque. We now discuss the implications 
of our results for the millisecond radio pulsars and for the 
possibility of gravitational waves from wobbling NS's. 

Ever since the discovery of the first millisecond radio 
pulsar (Backer et al. 1982), observers have been looking 
for rapidly rotating neutron stars near the breakup limit 
of v h w 1476 Hz M 1 1/ 4 2 J?~ 3/2 (see Cook et al. 1994 for an 
exhaustive survey of the breakup limits for different nuclear 
equations of state). However, few have been found. Even 
today, it is still the case that the majority of millisecond 
radio pulsars are spinning at 300 Hz or less. By consider- 
ing the characteristic ages of those pulsars that reside in 
binaries, Backer (1998) concluded that most of the initial 
spin frequencies of the radio pulsars are near 300 Hz, consis- 
tent with the later inferences from the accreting population. 
In this paper, we have shown that accretion-induced lateral 
density variations in neutron star crusts can indeed account 
for limiting the initial spins of millisecond pulsars. 

However, there are a few millisecond pulsars that 
do spin much faster than 300 Hz (e.g., B1937+21 and 
B1957+20, which have periods of roughly 1.6 ms). For these 
fast pulsars, a possible explanation is that their crusts are 
very uniform, so that the composition asymmetry is <C 10% 
as required to generate the appropriate Q22- A much more 
likely hypothesis is that these pulsars were spun up in tran- 
siently accreting systems with recurrence times greater than 
a few years, or accreting at a very slow rate. Lateral temper- 
ature variations in this regime would not persist, and hence 
the quadrupole moment resisting accretion would be much 
smaller. 

We finally turn to the possibility of gravitational waves 
from wobbling NS's. Our calculation of the maximum Q22 
sustainable by the crust can also be used to place a strong 



limit on the strength of such wobble radiation. A NS will 
"wobble" if its angular momentum J a is slightly misaligned 
from some principal axis of the crust. This phenomenon oc- 
curs in the Earth, and is called the Chandler wobble. It 
is a somewhat more complicated version of the torque-free 
precession of rigid bodies that is treated in undergradu- 
ate mechanics texts. For the sake of brevity, we will simply 
quote many of the relevant results from this area, and refer 
the reader to Munk & MacDonald (1960), Pines & Shaham 
(1972), Cutler & Jones (2000), and Jones (2000) for more 
details. 

Consider first a rotating (and hence oblate) NS, which 
is spinning along some principal axis of the crust. If the 
relaxed state of the star is spherical, then, no matter how 
the NS is kicked, it will not precess. However, if the crust 
has relaxed to the oblate state (i.e., the zero-stress state is 
oblate), then kicking the NS will cause precession. 

Let the star's angular momentum J a be along z a , which 
is displaced (by some kick) from the principal axis of the 
crust n a by angle 6 W . To a sufficient approximation, the 
inertia tensor of the wobbling body will be 

jab t ab . \ j / ~a ~b 1 ab\ . a j r ~a ~b 1 ab\ , * 7 * 7 \ 

I = J e +Ai(i(z z - -g ) + AI d (n n - -g ) . (77) 

and this "figure" will rotate rigidly around z a with preces- 
sion frequency ui p « fi. Thus AI d is the piece of the inertia 
tensor that "follows" the principal axis of the crust. Defining 
x a by h a « z a + 6 w x a , we see that I a b for the wobbling star 
can be re-written as 

r b = Ioe -\ {AIn+AId )^2 b -l g - b ) + w AI d (z a x b +x a z b ) .(78) 

Here the piece oriented along z a is "held up" by centrifugal 
force, while the piece proportional to 8 w AI d is "held up" by 
crustal shear stresses and is responsible for the star's Q21 
moment. A little algebra shows that 

Q21 = \f^e w AI d . (79) 

The dominant gravitational wave emission is at the wobble 
frequency, so tj gw w f2, and the gravitational wave luminos- 
ity is 

gobble = 2 ( ^ A/d) 2 n 6 = ^faSl* _ (80) 

Comparing with Eq. (3), we see that for a fixed spin fre- 
quency Q, 4= m=2 /L™ bble = 64(Q 22 /Q2i) 2 . The same ra- 
tio holds for the backreaction torques. Thus, if the accre- 
tion torque is balanced by the wobble radiation, rather than 
I = m = 2 radiation, resulting in v s ~ 300 Hz, then the re- 
quired quadrupole moment is 8 times larger for the wobble 
case. This would require an average strain a « 2.3 x 10~ 2 for 
M = 10~ 9 M Q /yr and our fiducial NS mass and radius. Al- 
ternatively, if one fixes the gravitational wave frequency (in- 
stead of the spin frequency), and asks, what is the maximum 
gravitational wave luminosity achievable by either I = m = 2 
radiation or wobble radiation, we see that the limit set by 
crust cracking is exactly the same for the two cases (up to 
small corrections of order {0,/Q.max) 2 and Q / M R 2 )-& rather 
pretty result! 
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Pressure (erg cm 3 ) 

Figure Al. Top panel: The run of the shear modulus fi/p with 
depth. Bottom panel: The run of T p (solid line) and |T Mo | (dashed 
line) with pressure in our model. Note that r Me is negative, so its 
absolute value is plotted. For pressures p < 10 30 erg cm -3 (i.e. 
for densities less than neutron drip) T p = |r^ e |. 



APPENDIX A: COEFFICIENTS AND SOURCE 
TERMS OF THE ELASTIC PERTURBATION 
EQUATIONS 

In this section we describe how to compute the coefficients 
and source terms of Eqs. (43) from the background models 
developed in § 2 and 3. 

We use the shear modulus as computed by Strohmayer 
et al. (1991) by Monte-Carlo simulations of both bec crys- 
tals and quenched solids. Their results can be conveniently 
rewritten in terms of the pressure of degenerate relativistic 
electrons, 



JL - 6 x 1(T 3 (Z\ 

p~ e ~ l + 0.595(173/r C oui) 2 V 8J 



2/3 



(Al) 



where Tcoui is defined in Eq. (19). We neglect the slight 
dependence of the shear modulus on temperature. Because 
of the uncertainty of the charge of nuclei at the bottom of 
the crust, we varied the numerator of Eq. (Al), as described 
in § 5.1. At densities higher than neutron drip, we have 
ai = (jj,/p e )(p e /p). The run of [i with pressure in the crust 
is shown in the top panel of Figure Al. 

In addition to the shear modulus, we require a few ther- 
modynamic derivatives. Let us first consider the case of the 
perturbations generated by a smooth composition variation 
in the crust, i.e., /j, e (r,6,<f>) = £t e (r) + A/j, e (r)Yi m (6 , (f>). In 
our computations we took Afj, e (r)/ fj, e (r) = const, but the 
formalism applies for arbitrary variation of Afi c with depth. 

In this case, we write p — p(p,fj, e ) and Ap using 



Eq. (34). The index F in this case is computed at constant 
composition, 



d hip 



<91np 



4Pe 

3 p + 



dlnp n 



dlnp 



: 

P 



(A2) 



where p e oc (p/p ) 4//3 is the relativistic electron degener- 
acy pressure, and p n (p,X n ) — p n {X n p/mb) is the neutron 
pressure. If the neutrons were completely degenerate, non- 
relativistic and non-interacting, the prefactor of p n /p would 
be 5/3. The actual run of Y p with pressure is shown by the 
solid line in Figure Al. 

Following Eq. (45), we write the source term as AS* = 
r Me (A^ e //i e ), where 



dlnp 



d In fi c 



3? 



cHnp„ 



d In /Lie 



P 



(A3) 



Using the chain rule, we express the prefactor of the p n /p 
term in terms of the prefactor of the similar term in Y p as 



d\np n 



d In fi c 

where 
dlnX n 
din fi c 



d\nX n dlnp n 



din n e dlnp 



X„ 



(A4) 



X„ 



Xnlle 

i-x n 
x a 



■ X n l 



Zi 

A! 



Z2 

A 2 



(f 



X n l 



(A5) 



which follows from Eqs. (14) and (15). Using the above for- 
malism, one can compute the source term AS and the index 
r of Eq. (45). The source term AS for this case is shown in 
the bottom panel of Figure 13. 

Most of the discussion in this paper was concerned with 
the quadrupole moments generated by shifts in capture lay- 
ers. For most of the calculation we have adopted the physical 
picture described by Eq. (46), i.e., that a temperature varia- 
tion 5T "undoes" some electron captures on the "cold" side 
of the star and induces more captures on the "hot" side, 
and hence moves the capture layer vertically. This picture 
is physically consistent, and the formalism for deriving the 
source terms is described below. However, it is somewhat 
easier to first understand the Lagrangian A« d picture intro- 
duced in § 1.2. 

In this picture, we first imagine that the crust is in- 
finitely rigid, and hence does not allow either vertical or 
horizontal elastic adjustment. We then go in and "by hand" 
undo some captures on one side of the star, to a height Azd 
above the original capture layer, and induce some captures 
to the same depth on the other side of the star. The per- 
turbed crust then looks like the sketch in Figure 2. After 
this perturbation, we allow the crust to respond elastically. 
This perturbation is Lagrangian in the sense that we know 
a priori which fluid elements have had their fj, e perturbed. 

We now formulate the above picture mathematically. 
The pressure at which the capture reaction occurs varies 
laterally. In a spherically symmetric star, the capture layer 
occurs at r — r c (i.e., (Ai, Z\) element is, say, 50% depleted 
at r = r c ). When we put in the fj, e perturbation, the deple- 
tion occurs at r — r c + Az^Yi m . It is an excellent approxima- 
tion to assume that the shape of the capture layer does not 
change, only its position in the star. In this case, we write 
p c (r, 9, 0) w fi c (r - r c - Az d Yi m (9, <f))) and Taylor-expand 
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to find A/i //i « — (d\n fi c /dr)Az d Yi m . The source terms 
arc then 



AS : 



-IV. ^A*. 



(A6) 



Mo dr 

where /ti e (^) is the run of electron mean molecular weight in 
the equilibrium model. Since p e (r) is constant outside the 
capture layers, this source term is non-zero only in the cap- 
ture layers. A plot of din p, c /dlnr is shown in Figure A2. 
Since He{r) resembles a step function with steps at the lo- 
cations of the capture layers, dln/i c /dlnr resembles a delta 
function at the locations of the capture layers. 

Now consider the case where the capture layers are 
moved around by lateral temperature variations, i.e. the case 
of 5T-sourced perturbations. In our equation of state, the 
dependence on temperature only comes in through p c , i.e. 
p = p(p,fi e {p,T)). Therefore, 



dlnp 



d Inp 



dlnp 



d In /i c 



dlnp c 



dlnp 



(A7) 



The derivative at constant p in the second term above is 
displayed more conveniently by changing back to (p, p c ) as 
independent variables, 



dlnp 



din p c 

so in case of ST perturbations, 
dlnp c 



r = r 



i-r B 



dlnp 



(A8) 



(A9) 



The term in square brackets is equal to 1 outside the capture 
layers, and hence F reduces to the same value as in the 
case of the perturbations caused by a smooth composition 
gradient (see Figure Al). However, the term in the square 
brackets, which we compute below, is significantly smaller 
than 1 in the capture layers. 

Using the chain rule as in Eq. (A7) and using Eq. (A8), 
we find that 



dlnp 



dlnT 



r 



d In fi c 



dlnT 



(A10) 



We can now evaluate the prefactor of the source term by 
switching independent variables from (p, T) to (p, T) to ob- 
tain 



dlnp 



dlnT 



din p c 



dlnT 



i-r„ 



d In p c 



dlnp 



(All) 



Eqs. (A9) and (All) allow us to evaluate both the source 
term AS and the index F (see Eq. [46]), if we can find the 
partial derivatives of p e (p,T) with respect to p and T. We 
now describe how to find these derivatives. 

While constructing an equilibrium stellar model, we 
computed the dependence of p c on r for each capture layer. 
However, at least in principle, we could compute p e (p,T) 
using the following approach. We combine the rate equation 
(17) and the equation of hydrostatic balance (18) to obtain 

dlnX-i 



dlnp 



V 
mg 



R ec (p,X, T), 



(A12) 



where rh — M /Aitr 2 is the local accretion rate. Approximat- 
ing r 2 and g as constant over the capture region, we then 
solve Eq. (A12) to obtain Xi(p,T), which, from Eq. (15), 



immediately gives us" p c (p,T). Eq. (A12) was derived as- 
suming the neutron star is in hydrostatic balance, but we 
shall use it to determine p, e (p, T) even in the perturbed star. 
This is consistent, because corrections to (J, e (p, T) are due to 
the nonspherical perturbation itself, and so would enter our 
perturbation Eqs. (43) only as higher-order corrections. 

In practice, rather than to first compute fj, e {p,T) and 
then differentiate it, it is much more convenient to use a 
simple approximation to evaluate the derivatives of /j, e . This 
approximation is essentially the same as the one used to 
evaluate the source term in the case of Azd perturbations, 
Eq. (A6). We presume that the shape of the capture layer 
does not change as we perturb the local temperature by ST, 
but only its location shifts. Mathematically speaking, we 
write p e {p, T) « p c (p—p c (T)), where p c is the (temperature- 
dependent) pressure at the "center" of the capture region, 
which we define more precisely below. We use the tilde on 
jl c to distinguish this "approximate version" from the actual 
function /i c . 

With this approximation, 



din /i c 



dlnp 



din jl c 
dlnp 



1 din p c 
V dlnr 



(A13) 



where p c (r) is just the run of p a in the background model 
(see Figure A2). V is defined in Eq. (44) and is equal to 
r/h, where h is the local pressure scale height. In the crust 
h/r <J 0.1 (<C 0.1 in the outer crust). It is then evident from 
Figure A2 that the term in square brackets in Eq. (A9) is 
indeed much smaller than 1 inside capture layers. 
Similarly, we have 



d In p c 



dlnT 



dlnp, c dlnp c _ 1 dlnp e dlnp c 



(A14) 



dlnp dlnT V dlnr dlnT 

In the above equation, p c is a function of T for a given 
capture region, and is not a function r or p. 

Eq. (A12) exhibits the competition between the local 
compression timescale t comp = p/rhg and the electron cap- 
ture timescale t cc = 1/R CC . Most electron captures happen 



when t c 



t e 



We use this to define the center of the 



capture layer, p c : 

^R ec (E F (p c ),T) 
mg 



(A15) 



Differentiating the above equation with respect to T and 
using the capture rate (11) we find 



din p c 
dlnT 



E F -Q 



k B T 



1 + 



E F dlnEp 



k B T dlnp 



(A16) 



All the terms on the right-hand side of the above equation 
are evaluated at p = p c . But to find p c we would have 
to solve the transcendental Eq. (A15). However, we note 
that the above expression enters in Eq. (A14) multiplied by 



S3 To be more precise, fi c is a function of p and a functional of 
T, since it depends not just on the T at the point but on T in 
the whole capture region. However, Bildsten & Cumming (1998) 
showed, and we confirm for the case of deep capture layers that 
include free neutrons, that the variation in T over the capture 
region is quite small. Thus we treat T(r) as constant over the 
capture region, in which case there is no practical distinction be- 
tween a function and a functional. 
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Figure A2. Left (a): The run of the mass fraction X\ (top panel) and din fi e /dlnr with radius for a shallow capture layer at threshold 
energy Q = 23 MeV. Right (b): Same as (a), but for a deep capture layer at Q = 95 MeV. 



rfln/i e /dlnr, which is a sharply peaked function of r. Thus 
we can treat the right hand side of (A16) as a function of 
r and let the delta-function like shape of dln/i e /c?lnr pick 
out the correct value of p ~ p c - 

Finally, to evaluate Eq. (A16) we need 



dln-Ep _ 1 
d\np 3 



1 1 



+ 1 



din n c 



dlnr 



(A17) 



which follows from Ef oc (p//^) 1 / 3 and a modest amount 
of algebra. We can now evaluate both F and the source 
term AS using the properties of the background model: use 
Eqs. (A9) and (A13) to compute F, and Eqs. (46), (All), 
(A13), (A14), and (A16), as well as the local 5T, to compute 
AS. The source term AS is plotted in the bottom panels of 
Figures 10-12. 

Using Eq. (A16) we can establish a precise relation be- 
tween the vertical shift Aza of the "center" of the e~ capture 
layer (in absence of elastic readjustment) and the local tem- 
perature perturbation 5T, 

Az d 1 dln Pc ST 

— = ^^hTr—- (A18) 

This relation is plotted in the bottom panel of Figure 3 in 
§ 2.1. Similarly, the change in electron Fermi energy at which 
most captures occur in response to a local temperature per- 
turbation is just 

dlnp c dlni? F E F 



AEf 



dlnT dlnp k B T 



k B 8T, 



(A19) 



or AEf = -r(k B 8T), where 



T = 



d In pc d In Ef Ef 



(A20) 



llnT dlnp k B T' 

The function T is plotted in the top panel of Figure 3 in 
§ 2.1. The complicated dependence of Azd and T on the 
capture layer depth stems from the differences in the number 
of electrons captured and number of neutrons emitted in 
each capture layer. 
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Table Al. List of Variables 

Name Description Section 

M Mass accretion rate 1 

u s Spin frequency 1 

B Magnetic field 1 

Q22 Mass quadrupole moment 1 

-^Edd Eddington accretion rate 1 

N a Accretion torque 1.1 

Sp Eulerian density perturbation 1.1 

Ql m Mass multipolc moment 1.1 

v gw Gravitational wave frequency 1.1 

N gw Gravitational wave torque 1.1 

M NS mass 1.1 

R NS radius 1.1 

Q eq Mass quadrupole moment needed for 1.1 

equilibrium between accretion and GW 

emission 

Mi. 4 NS mass in units of 1.4 M 1.1 

-R6 NS radius in units of 10 km 1.1 

I a b Mass quadrupole tensor 1.1 

Q Spin frequency 1.1 

a Average shear strain 1.1 

fs,eq Equilibrium spin frequency 1.1 

d Source distance 1.1 

E gw Energy loss due to GWs 1.1 

h a GW strain 1.1 

F x X-ray flux 1.1 

A Nuclear mass number 1.2 

Z Nuclear charge 1.2 

ST Lateral temperature perturbation 1.2 

Azd Vertical shift of a capture layer 1.2 

Ap Lagrangian density perturbation 1.2 

T Sensitivity of e~ capture layer location to 1.2 
temperature 

p Shear modulus 1.2 

&b Boltzmann constant 1.2 

p nt j Neutron drip density 2 

K Conductivity 2 

icomp Compression timescale 2.1 

fn Local accretion rate 2.1 

g Local gravitational acceleration 2.1 

e Electron charge 2.1 

rrtt, Baryon mass 2.1 

E-p Electron Fermi energy 2.1 

Q e~ capture threshold energy 2.1 

ft ft value 2.1 

_R oc e _ capture rate 2.1 

h Scale height 2.1 

t cc e~ capture timescale 2.1 

Xi Mass fraction of the e~ capture element 2.2 

X2 Mass fraction of the e - capture product 2.2 

X n Mass fraction of free neutrons 2.2 

n e Density of electrons 2.2 

p c Electron mean molecular weight 2.2 

v Compression flow speed 2.2 

AEp Thickness of a capture layer 2.2 

r Radius 2.2 

p e Electron pressure 2.3 

pn Density in units of 10 11 gem" 3 2.3 

n n Free neutron density 2.3 

M r Mass enclosed within radius r 2.3 

a Internuclear spacing 2.3 

Tcoul Coulomb energy density relative to k^,T 2.3 



Table Al - continued 
Name Description Section 

F Heat flux 2.4 

F r Radial component of the heat flux 2.4 

e nuc Energy generation rate 2.4 

e„ Urea neutrino emissivity 2.4 

e Energy gain/loss rate per gram 2.4 

_B nuc Total energy released in a capture layer 2.4 

A Superfluid gap energy 2.4 

pnuc Nuclear density 2.4 

Lcorc Energy lost from the NS core 2.4 

Tj-mm H/Hc burning temperature 2.4 

/ nuc Asymmetry in the nuclear energy release 3.1 

SE nuc Difference total energy released in a cap- 3.1 

ture layer 

rife Temperature sensitivity of the 3.1 

conductivity 

n e Temperature sensitivity of the neutrino 3.1 

emissivity 

/comp Asymmetry in conductivity 3.1 

S(Z 2 /A) Difference in charge-to-mass ratio 3.1 

SF a Flux perturbation 3.2 

f a unit radial vector 3.2 

V a Derivative operator 3.2 

K COTe Conductivity of the core 3.2 

_ftr crus t Conductivity of the crust 3.2 

SF T , core Radial flux perturbation in the core 3.2 

SF± coro Transverse flux perturbation in the core 3.2 

Ail Crust thickness 3.2 

Fa Accretion flux 3.3 

/ Composition or nuclear heating 3.3 

perturbation 

t t h Thermal time 3.3 

£ a Displacement vector 4 

Qf^ 1 "Historical" multipole moment 4 

Qf^ Multipolc moment due to the "current" 4 

perturbation 

Ap e Smooth composition gradient 4 

r ao Stress tensor 4.1 

g a \y Metric tensor 4.1 

£ r Radial component of the displacement 4.1 

£x Transverse component of the displacement 4.1 

$ Gravitational potential 4.1 

(i + 1) 4.1 

T rr , T r x Components of the stress tensor 4.1 

e a b, fab, Aa6 Angular tensors 4.1 

zi — Z4 Eigcnfunctions of the elastic perturbation 4.1 

equations 

ai — 04 Coefficients of the elastic perturbation 4.1 

equations 

r Compressibility 4.1 

AS Source term 4.1 

U, V Polytropc variables 4.1 

ftop, »*bot Top and bottom of the crust 4.2 

p s ,pi Density of the solid and of the liquid at the 4.2 

crust-core boundary 

Q22 Plane-parallel approximation to the 4.2 

quadupole moment 

a + max Breaking strain of the crust 6 

t a b Stress tensor (nonlinear) 6.2 

t rr , t r ±, tjy Components of t a b 6.2 
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Table Al - continued 



^ab 


O Lr! <X1±1 UL^llOLJJ. 


6.3 


(j 


Average strain 


6.3 


cy . . , rr .. i «t a 


Components of lTq,^ 


6.3 


ro 


Pressure at the bottom of the crust 


6.3 


Pb 


l~)pVI<^itv Pit flip bnttOTTI of tVlP P"T1T>t 

Li V QiL Li±l\_. UUljlrVJlll Ul trills U.O Lr 


6.3 


(u/v) 


A vpT"a o"p ^Vi pa r rnnrl n 1 ii^ 


6.3 




Effective strain 


6.3 


crust 


Moment of inertia of the crust 


6.5 




IVTritn PV1T PIT IVlPftlPl PIT ("MP TIPIirTPlTl GTPf 


6.5 


5<E> 


Gravitational potential perturbation 


6.5 




(^nvrppfinri Hup to tlip npclppt of <>plf- 


6.5 




pravitv 




J° 


Angular momentum vector 


7 




Direction of the rotation axis 


7 


n° 


Pricipal axis direction 


7 




Wobble angle 


7 


jab 


NS inertia tensor 


7 


I 


Spherically symmetric part of 


7 


Aln 


Oblate part of I a ^, for crust whose un- 


7 




stressed state is spherical 






Oblate part of I , for a fully relaxed crust 


7 


V21 


in — 1 mass Cjuadrupolc moment 


7 




Precession frequency 


7 




GW angular frequency 


7 


r wobble 


GW luminosity due to the wobble 


7 


Pc(T) 


Pressure at the center of the capture layer 


A 


fie 


Run of fj, e with radius 


A 




Logarithmic index of pressure with density 


A 




Logarithmic index of pressure with elec- 


A 




tron mean molecular weight 
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